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Abstract 



The present thesis provides algorithms for some problems in computational geometry 
using prune and search technique. We solve the intersection radius problems in the 
plane for sets of line segments, rays, wedges, half planes and convex polygonal 
disks. Then we compute a centerpoint of a planar set of points. This is followed by 
description of a generalised technique of prune and search. 

The intersection radius of a set of n geometrical objects is the radius of the smallest 
closed ball that intersects all the objects of the set. We have designed algorithms to 
optimally solve intersection radius problems for various kinds of objects. We show 
how the prune and search technique, coupled with the strategy of replacing a ray by 
a point or a line can be used to solve, in linear time, the intersection radius problem 
for a finite set of line segments in the plane. 

Next, the scope of this technique is enlarged by showing that it can also be used 
to find the intersection radius of a set of convex polygons in linear time. Moreover, 
it is immaterial if the set also contains other types of geometric objects like points, 
lines, rays, line segments, half planes and wedges. In fact, it is shown how such a 
mixed set of objects can be handled in a unified way; and this is the other important 
contribution of the thesis. Previously there existed no known algorithms to solve 
these intersection radius problems efficiently. 

The center of a set P of n points in plane is defined as the set of points c such 
that any closed half plane containing c contains at least |"n/3] points of P. A 
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centerpoint is any point in the center. It can be viewed as a generalisation of the 
median of a set of reals. In the thesis, it is shown how a centerpoint of a finite 
set of points in the plane is computed in linear time. The described algorithm is 
optimal which significantly improves the O(n log 3 n) complexity of the previously 
best known algorithm. We need to use suitable modifications of the ham-sandwich 
cut algorithm and the prune and search technique to achieve this improvement. 

The optimal sorting network by Ajtai et al. has a noteable feature that it approxi- 
mately sorts the input data in its intermediate steps. A technique is presented for 
solving computational geometry problems by exploiting this fact. This technique 
is a sequel to the Megiddo's technique of designing serial algorithms by applying 
parallel computation algorithms. We get optimal linear time algorithms for some 
problems by applying this technique. We do this by the synthesis of prune and 
search and parametric searching. 



Once there lived a village of creatures along the bottom of a great crystal 
river. Each creature in its own manner clung tightly to the twigs and rocks of 
the river bottom, for clinging was their way of life, and resisting the current 
what each had learned from birth. But one creature said at last, "I trust that 
the current knows where it is going. I shall let go, and let it take me where it 
will. Clinging, / shall die of boredom." 

The other creatures laughed and said, "Fool/ Let go, and that current you 
worship will throw you tumbled and smashed across the rods, and you will 
die quicker than boredom!" 

But the one heeded them not, and taking a breath did let go, and at once was 
tumbled and smashed by the current across the rocks. Yet, in time, as the 
creature refused to cling again, the current lifted him free from the bottom, 
and he was bruised and hurt no more. 

And the creatures downstream, to whom he was a stranger, cried, "See a 
miracle! A creature like ourselves, yet he flies/ See the Messiah, come to save 
us all/" And the one carried in the current said, "J am no more Messiah than 
you. The river delight to lift us free, if only we dare let go. Our true work is 
this voyage, this adventure." 

But they cried the more, "Saviour/" all the while clinging to the rods, making 


legends of a Saviour. 



Preface 



Geometry is one of the most ancient branches of mathematics. One of the earliest 
attempt to systematise Geometry was made as early as in 300 B.C. by Euclid in 
Greece. This field also has distinction of being the first ever subject to profess 
of being rational and logical. At least it was believed by the people at that time 
and even now of standing on strictly formal grounds of logic. The whole system 
in Elements by Euclid is built on five basic geometrical postulates and numerous 
theorems are proved using only these few postulates. It can be said without any 
doubt that Geometry is one of the oldest field and scores of brilliant mathematicians 
and researches have worked in this field enriching geometry considerably by their 
discoveries and findings. 

Computational Geometry on the other hand is only of recent origin. It deals with the 
computational aspects, its theory and applications. The subject matter of Compu- 
tational Geometry ranges from theoretical computer science to that of more applied 
nature in Algorithms. To categorise the field might be difficult but it can be safely 
presumed that Computational Geometry is an offshoot of Applied Mathematics. It 
profusely uses the results of Combinatorial Geometry and constructs of Euclidean 
Geometry in its theoretical framework and the results of Analytical Geometry in 
its applications. When it deals with proofs of existences, theorems concerning 
Geometry, classifications of geometrical objects in different equivalent classes or 
enumerations of these then we are concerned with theoretical Computational Ge- 
ometry and when it deals with computing geometrical objects, or their attributes 

vi 



or their count then we are concerned with application of Computational Geometry. 

It is important to know the types of problems in Computional Geometry to under- 
stand it better. The fundamental problems are that of construction of convex hulls, 
location of points in subdivisions, construction of Voronoi diagrams, etc. For a more 
complete list the pioneering doctoral work of Shamos can be referred. 

In this thesis we are concerned with the algorithmic aspects of Computational 
Geometry. We give algorithms to compute geometrical concepts such as intersection 
radius and center points. Intersection radius has got reference in the first ever 
collection of problems in Computational Geometry by Shamos where as the concept 
of center is dealt with in the book on geometry by Yaglom and Boltyanskii. The 
last part of the thesis deals with a generalised computational tool that can be used 
for geometrical optimisation problems. 

I wish to acknowledge my thesis supervisor Prof. Ashish Mukhopadhyay for his 
guidance during my Ph. D. work. I would like to acknowledge the work of researchers 
(especially Herbert Edelsbrunner and Nimrod Megiddo) before me in the field of 
computational geometry for providing motivation to persue this work. I thank my 
parents for leaving me alone to persue higher studies. I thank all my friends in and 
out of I.I.T. Kanpur for the company. I thank the I.I.T. Kanpur administration for 
bearing me. I thank the workers' cooperative of I.I.T. Kanpur for allowing me a 
chance to do social work. I finally thank myself for completing this thesis. 
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Chapter 1 



Introduction 



This thesis deals with a method of computing intersection radii for various kinds of 
geometrical objects and also of computing a centrepoint for a planar set of points. 
It also describes a general technique for applying the prune and search paradigm for 
solving geometric optimisation problems. 



1.1 Intersection Radius 

Suppose there are several sites in a two-dimensional metric space and we have to 
choose a point p such that it is "nearest" to all the sites in question. Here the 
word "nearest" implies that the point p minimises the maximum distance of the 
sites from p. If a circle is drawn to "span" (cover) all the sites then it is called a 
spanning circle. When we draw a spanning circle centered on the "nearest" point p 
with radius equal to its distance from the furthest site(s), then this circle is called 
the minimum spanning circle. 

Phrases such as "spanning circle" and "minimum spanning circle" are meaningful so 
long as we are interested in spanning only point sites. These phrases are inadequate 



to describe the problem of computing a smallest closed disk that intersects other 
objects such as lines, rays, line segments, etc. Thus the notion of stabbing is 
introduced. 

When one object intersects or touches another it is said that the former stabs the 
latter. An object is called a slabber of a given set of objects if it stabs each object 
of this set. For example, the closure of the interior of the minimum spanning circle 
above is the smallest closed disk that stabs the given set of points. 

Let C be a finite collection of geometrical objects in a d-dimensional Euclidean space. 
The stabbing problem consists of finding an object (the stabber) which intersects each 
member of C. Typically, the stabber could be a line, a hyperplane or a disk, and C 
could be a collection of points, lines, line segments, rays, hyperspheres, polyhedron 
or any mix of these. A survey of some recent results is available in the paper by 
Houleetal. [15]. 

A classical problrm in this area is that of finding the intersection radius of a finite 
set of points in the plane, which is also known as the 1-ccntre problem [30, 31, 6]. 
It was shown by Mcgiddo, and subsequently by Dyc-r [22, 12, 11], how this can be 
solved in linear time for fixed d. However, until a recent paper [3] no attempt was 
made to extend this to a more complex set of objects than points or to a collection 
containing different kinds of objects. 

Attempts have been also made recently to find more complicated slabbers for the 
stabbing problem, or to find the "best" stabber which optimises some measure 
defined on the class of stabbers in question. Goodrich and Snoeyink [14] presented 
an 0(r?logr?) algorithm to find a convex polygon whose boundary intersects each 
of T? parallel line segments. Rappaport and Meijer [24] showed that a perimeter 
minimising polygonal disk that intersects each of n parallel line segments can be 
found in O(nlogrc) time. They have also extended their result to a set of isothetic 
line segments. Mukhopadhyay and Kumar [27] have shown that for a set of parallel 
line segments an area minimising polygonal disk can also be found in O(nlogn) 
time. Bhattacharya and Toussaint[4] gave an 0(nlog 2 n) algorithm for computing 



the shortest line segment that intersects a set of n given line segments in the plane. 
This bound was later improved to O(nlogn) by Bhattacharya et al.[2]. 

In the present thesis, the stabber is a disk. The stabbing problem in this case is 
known as the intersection radius problem[13]. The intersection radius of a planar 
collection of objects is the radius of the minimum stabbing disk of this collection. 
There is no known algorithm which solves the intersection radius problems for line 
segments. There are some variants of this problem which have been the subject 
of research of various researchers. As for example, when the stabber is a vertical 
line segment and the objects to be intersected are lines, the intersection radius 
problem is simply the Chebyshev approximation of points in the dual space[5, 29). 
This problem can easily be solved in linear time by transforming it into a linear 
programming problem. 

In this thesis, the intersection radius of a collection of line segments is computed 
by combining the prune-and-search strategy of Megiddo [21] with the strategy of 
replacing line segments with points or lines [3]. The scope of this technique is 
enlarged by showing that the intersection radius of a collection of convex polygons 
can also be computed in linear time. Further, it is immaterial if the collection 
contains a mix of other geometric objects such as lines, points, rays etc. for we show 
that it is possible to treat such a mixed collection of objects in a unified way. 



1.2 Centrepoint of a Planar Set of Points. 

We all have an intuitive idea as to what phrases like "the very center of the square" 
or "the very center of the city" mean. To capture this intuition in a quantitative 
way, the center of a set of n points, P, in 9^ is defined as the maximal subset of W* 
such that any closed halfspace intersecting this subset contains at least \n/(d -f 1)] 
points of P [35]. This subset is non-empty for any finite configuration of points (see, 
for example, [13]). Furthermore, it is closed and convex. A centrepoint is a member 
of the center of P. 



On the real line 3?, a centrepoint is no other than a median of P. Thus a centrepoint 
can be viewed as a generalisation of the median of a set of reals. On the other hand, 
the center can also be viewed as a particular fc-hull of P. The k-hulloi P is a maximal 
subset (closed and convex) of $t d such that any closed halfspace intersecting this 
subset contains at least k points of P. For instance, the 1-hull of P is its convex hull 
and the center is its \nf(d+ l)~|-hull. The property of balanced partitioning makes 
the centrepoint useful for efficient divide and conquer algorithms in geometrical 
computing and large scale scientific computing[26, 32, 25, 36]. Recently Donoho and 
Gasko have suggested that centrepoint can be used as "robust" and high "breakdown 
point" estimators for multivariate datasets [10]. 

The interesting algorithmic problem of computing a centrepoint has been considered 
by various researchers. Cole et al. gave an 0(nlog 5 n) algorithm for computing a 
centrepoint of a planar set of points [9]. Subsequently, Cole improved this bound to 
0(nlog 3 n), using the powerful technique of slowing down a sorting network [8]. In 
this paper, we propose an optimal linear time algorithm for computing a centrepoint 
of a planar set of points by using suitable modifications of the ham-sandwich cut 
algorithm for a pair of separable point sets [23] and the prune and search technique 
of MegiddoplJ. 

Linear time algorithms, however, were known for computing an approximate or e- 
centrepoint [19, 32, 23]. We obtain this weaker type of centrepoint if we decrease 
the lower bound, in the above definition of the center, to |~7?(1 e)/(d -f 1)1, where 
< < 1. Actually, Megiddo [23] only gave an algorithm for computing a partition 
of a (planar) set of n points with two lines such that each closed quadrant contains at 
least [n/4j points. An algorithm for computing an e-centrepoint, where < < 1/4, 
is implicit in this. 

The thesis proposes an optimal algorithm for computing a centrepoint of a planar 
set of points by using an interesting modification of Megiddo's prune-and-search 
technique[21). This consists of adding a few extra points in each pruning step so that 
a subspace of the original solution space is retained, while ensuring a net deletion 



of points. In the description of our algorithm, we assume the usual RAM model of 
computation; the point set P, however, is not assumed to be in general position. 



1.3 Prune and search in Slowed Down Sorting 
Networks 



A substantial part of Computational Geometry deals with designing efficient algo- 
rithms for given problems. More often then not, the solutions of these depend on 
the ingenuity of the inventor of the algorithm. However, there are also attempts 
to provide general tools to solve the problems of a specific type. Paradigms like 
divide and conquer, parametric searching, line sweep, prune and search, linear 
programming etc. go a substantial step towards this goal. 

The present thesis provides such a tool. We modify the method of parametric 
searching and slowed down sorting networks by Cole [8] to get a new technique. 

The parametric searching of Megiddo [20] is as follows Let there be an efficient 
parallel algorithm for a problem A such that solution of A can be used in the solution 
of another problem B. Then, in some cases, we get an efficient sequential algorithm 
for B by exploiting the efficient parallel mechanism of the parallel algorithm for 
A. This technique has been applied to a wide variety of problems yielding efficient 
algorithms. In particular, we achieve good results for the parametrised problems 
that use parallel sorting algorithms in their solutions. We do this by replacing 
the evaluation of the parallel comparisons of an iteration in the parallel version 
by simultaneous resolution of these in the serial version. The running time of these 
algorithms is further improved by the introduction of weights in the comparisons [8]. 
We then simultaneously evaluate at least a fraction of total weight of these in every 
iteration to design more efficient algorithms. 



We further improve the running time of the above algorithms by introducing, wher- 
ever applicable, prune and search in these. We do this by seeking to compute 
the fc-th largest element of the input set instead of seeking to sort it as in the 
previous techniques. For this, we run the sorting algorithm on a given input for a 
few iterations and then prune the set. It can be easily seen that this approach is 
useful only where pruning is applicable. 



1.4 Organisation of Thesis 

The organisation of thesis is as follows. In chapter 2 we solve the intersection radius 
problem for planar line segments in linear time. This algorithm is applied to all 
kinds of objects such as wedges, rays, halfplanes and convex polygonal disks in the 
next chapter. Moreover, it is shown how such a mixed set can be treated uniformly. 
In chapter 4, we present a linear time algorithm for computing a centrepoint of a 
finite planar set of points. We present a generalised tool for applying prune and 
search to a given problem in chapter 5. We conclude and discuss directions for 
further research in the last chapter. 



Chapter 2 



Intersection Radius of a Set of 
Line Segments 



The intersection radius of a finite set of geometrical objects in a ^-dimensional 
Euclidean space, E d , is the radius of the smallest closed ball that intersects all 
the objects of the set. In this chapter, we describe optimal algorithms for some 
intersection radius problems in the plane. We first briefly review a linear time 
algorithm by Megiddo to determine the smallest disk that intersects a finite set 
of points. Next we show how the prune and search technique, coupled with the 
strategy of replacing a ray by a point or a line can be used to solve, in linear time, 
the intersection radius problem for a finite set of line segments in the plane. 



Previously, no efficient algorithm was known to solve the intersection radius problem 
for line segments in the plane. 





() 



Figure 1 Spanning circles and stabbing disks of points and line segments 

2.1 Introduction 



Suppose there are several sites in a two-dimensional metric space and we have to 
choose a point p such that it is "nearest" to all the sites in question. Here the word 
"nearest" is to be* understood in the following way: the nearest point ;> minimises 
the maximum distance from the sites to p. This problem frequently arises in practice 
in the installation of centralised facilities e.g. such as installation of a transmission 
centre. If a circle is drawn centred at p with radius equal to the distance from p to 
the farthest site, to span all the sites, then it is called a minimum spanning circle 
of the sites (Fig. l(a)). 

The phrase "spanning circle" is meaningful only if we are interested in spanning or 
covering points and is inadequate to describe the problem of computing the smallest 
closed disk that intersects other objects such as lines, rays, line segments, etc. Thus 
we introduce the notion of stabbing for these objects. 

When one object intersects or touches another we say that the former stabs the 
later. A slabber of a set of objects is an object that stabs each member of this set. 
For example, the closure of the interior of the minimum spanning circle above is the 



smallest stabbing disk of the given set of points. 

The concepts of spanning circles and stabbing disks, though related, are different 
in several ways. To span a set we only need to consider extreme points of the 
individual members of the set. For example, to span a set of line segments we only 
need to compute a spanning circle for the cndpoints of these. However, this does 
not help in the computation of minimum stabbing disk of line segments which is 
frequently smaller. Similarly, in general also, it docs not help to characterise the 
minimum stabbing disk if the minimum spanning circle of the given set is computed 
(Fig. l(b)). 

The intersection radius of a planar set of objects is the radius of the minimum 
stabbing disk of this set. Now this question naturally arises how to compute the 
minimum stabbing disk/intersection radius for a given set of objects efficiently in 
linear time. This problem is progressively solved in this chapter and the next, first 
for simple objects, viz points and lines, and thru for more complex ones. 

The problem of computing the intersection radius usually reduces to a non-linear 
programming problem except when the stabbed set consists of straight lines alone. 
Not only for these but also for theii countc-i parts in higher dimensions, i.e. hyper- 
planes, the intersection radius problem similarly reduces to a linear programming 
problem [3]. For fixed d, this latter problem can be solved in linear time, using the 
algorithm of Dyer or Megiddo or Olarkson [11, 21, 7). Or instead, a recent simpler 
arid more efficient randomised algorithm of Seidel [28] can be used that runs in 
O(?7 (d + 1)!) expected time where n is the number of constraints in the linear 
piogramming problem. 

The case of points also has been fully tackled [21]. This was done by exploiting the 
fact that the bisectors of points are straight lines. These bisectors can be seen as 
linear functions and therefore the methods of linear programming problem had been 
applied to this case too. 

The organisation of this chapter is as follows. In Section 2.2 wo give a few definitions 
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arid discuss the motivation of the method to compute the intersection radius of points 
in linear time. The solution of the intersection radius problem for points in plane 
is presented in Section 2.3. In Section 2.4 we design algorithms to compute the 
intersection radius of lines and rays which are used to compute intersection radius 
of line segments in plane. 



2.2 Preliminaries 

Let P be a set of n points in the Euclidean plane and let c be the centre of a 
spanning circle of radius r. If d(c,p) denotes the distance of c from a point p in P, 
then clearly, 

r 



, 

Denoting by, d(c,P), the distance of c from the set P, which is also the radius of 
the smallest spanning circle centred at c, we have 

d(c,P) = max</(c,p) 

pP 

The intersection radius, IR(P), of P is the radius of the minimum spanning circle. 
Thus, 

1R(P) = mmd(c,P). 

We can similarly define the intersection radius of a set of lines, rays, line segments 
etc. We will denote minimum spanning circles by min-circle and minimum stabbing 
disks by mm-disk . The centres of these in the unconstrained and constrained 
versions of the problem will be denoted by unconstrained centre and constrained 
centre respectively or just centre when the context is clear. 

Let us first see how intersection radius of P is computed. Consider the min-circle 
of three points in the plane. The circumscribed circle of the triangle formed by 
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p 
Figure 2: Computation of intersection radius of T> 

these points is also a spanning circle of these. However, it is not necessary that 
this circle is a minimum too, unless the triangle formed by the points is an acute or 
right-angled triangle. In the case of obtuse-angled triangle the min-circle is obtained 
by drawing a circle with its diameter as the largest side of the triangle. A naive 
algorithm of complexity O(n 3 ) to compute the intersection radius of P now is same 
as computing the radius of the largest of rnin-circles of all subsets of P of size three. 
It can be shown that the largest of those min-circles spans P. 

A better algorithm, that runs in O(n 2 ) time, can be designed as an improvement to 
this naive algorithm by fixing two points in each iteration such that these two points 
are on the circumference of the spanning circles. These points are then dynamically 
updated until the min-ciiclo of P is finally obtained. The details arc- as follows. We 
start with two fixed points, p, q P, such that the points in P lie on the same side 
of pq. The min-circles of three points are drawn for p, q and every other point of 
P. Let the point corresponding to the largest of these circles, which incidently also 
spans P, be r. If Apqr is obtuse-angled then we delete the obtuse-angled vertex 
and repeat the same steps with the largest side as chord, otherwise the minimum 
spanning circle of A/x/r is the min-circle of P(Fig 2). 

We might think of improving this algorithm further by fixing only one point and 
then computing the spanning circles. But, since we need two more points to fix a 
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spanning circle, and we have a quadialic numbei ol choices, sonic other approach 
is needed. We therefore consider fixing a line on which the centres of the spanning 
circles are constrained to lie. Further, we also have to think of some way of not 
considering all of the quadratic choices that determine a spanning circle. This leads 
us to the important idea of prune and search. 

The radius of the min-circle of P such that centres are constrained to be any point 
p is a convex function of p, f(p). In particular, when the domain of this function 
is restricted to a line, /() is still convex. Therefore by examining the gradient of 
/() at any arbitrary point p on a query line we can tell whether f(p) is minimum 
and if /(;>) is not minimum then on which side of p on the line the minimum of 
/() is attained. We will see in the next section how these observations help us in 
computing intersection radius of points in linear time. 



2.3 Intersection Radius of Points in the Plane 

When the objects to be intersected are only points, the algorithm of Dyer or Megiddo 
[12, 22] for the 1-centre problem can be used. A brief account of the algorithm is 
presented in this section, such that, it can be further generalised to include other 
kinds of objects. 

Let P be the set of points whose intersection radius is to be determined and let its 
size be N p . 

To begin with the following question is answered: Which are the points that can be 
pruned. To answer this, we note that min-circle is determined by a set of only two 
01 three points. These points are farthest from the centre of min-circle among all 
the points in P. We try to determine these by removing the other points of P in an 
organised manner as follows: 

// the centre of min-circle is contained in a half plane determined by the 
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Figure 3- Localisation of Centre of Mm-disk of P 

bisector of some two points of P then th( point which is nearer to the 
centre can be dropped from P in the computation of intersection radius. 

To prune away a constant fraction of points according to this rule, we need to localise 
the centre of min-cirde with respect to a constant fraction of bisectors. The crux of 
the method is How can this be done? 

We further make note of the fact that, given a non optimum spanning circle of 
radius r, the radius of min-< iide is smaller than r and therefore, the centre of min- 
circle will be at a distance smaller than r from any of the points in P. This has 
the following consequences. First, for any point P, we can localise the centre in 
a small conical region whose apex is P. This can be done by first computing the 
smallest possible spanning circle centred at P and then computing intersection of 
interior of cirdes with radii equal to the radius of this circle centred at points on 
its circumference, of which at most three need to be considered for the purpose at 
hand (Fig. 3). Secondly, for any line L, we can localise the centre in one of the sides 
of Z/, by similarly computing the min-circle whose centre lies on L first and then 
computing the conical region, which can be shown not to intersect /,, with respect 
to the centre of this circle. 

To compute intersection radius of P in linear time we consider the following related 
subproblcms. Subproblcm 2 and Subproblem 3 can be solved using oracles of 
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Subproblem 1 and Subproblem 2 respectively. The convexity of the distance function 
is crucial in both the solutions as it also provides the direction, consequently a half 
space, in which the search for the centre of min-circle is localised. 

Subproblem 1 Compute the min-cjrcle ofP such that the centres of the spanning 
circles arc at a fixed point P. 

Subproblem 2 Compute the mi u -circle ofp such thai Iht ant res of the spanning 
circles lie on a fixed line L. 

Subproblem 3 Compute the min-circle of P. 

The solution to the first subproblem is the simplest. We compute the maximum of 
the distances of points in P to the fixed point P and choose this distance as the 
radius of min-ciicle. This computation takes linear time. 

The min-circle, whose centre is constrained to lie on line L, is determined as follows. 
We pair the points arbitrarily in [Ap/2j P a i rs anc * compute the perpendicular 
bisectors of each pair. If any bisector is parallel to the line L then we drop the 
point that is nearer to L in the corresponding pair. In case the bisector is same as 
line L then any one of the points in the pair is dropped. After this, Subproblem 1 
is solved for the median point, P, of the intersection points of L with the rest of 
bisectors. If /' is the constrained centre then the computation is over, otherwise 
one point corresponding to each of these bisectors can be dropped by localising the 
constrained centre on L with respect to P. Thus it can be seen that at least [/V p /4j 
points are dropped in one iteration. We iterate until fewer than four points remain 
when any straight forward algorithm is applied. This is the technique of prune and 
search in which the size of input set is truncated by a constant fraction in each 
iteration. It is easily seen that it runs in linear time in the present case. 

The solution of Subproblem 3 is obtained as follows. We pair the points in P and 
form their perpendicular bisectors. Let the median slope of non-vertical bisectors of 
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these be s m . The bisectors of slope strictly larger than s m are paired with bisectors 
of slope strictly smaller than s m . As it can be seen, this pairing may not cover all 
the bisectors, if the number of positive slope bisectors is not equal to the number 
of negative slope ones. Next, intersection points of each of these bisector pairs are 
computed. 



First we localise the centre of min-circle with respect to a vertical line which divides 
the total of the vertical bisectors and the intersection points above in two halves. If 
the centre of min-circle lies on this vertical line then the min-ciicle is determined, 
otherwise let J y be the halfplane, determined by this line, in which the centre of 
min-circle is localised. Next we localise the centre of min-circle with respect to a 
line with slope s m such that it divides the total set of the intersection points in J y 
and the bisectors with slope s m in the above. As previously, if the centre of min- 
circle lies on this line then we terminate our computation, otherwise let J m be the 
halfplane determined by this line in which the centre of min-circle is localised. Then 
we drop a point which is nearer to the centre of min-circle for each corresponding 
bisector that does not intersect J v H J m 



It can be easily shown that there are at least [W p /4j bisectors which are one 
of these the paired bisectors, the bisectors with median slope or the vertical 
bisectors. Thus, in every iteration at least |_/V p /16j points are dropped. These steps 
are repeated with the truncated sot of points until there aro fewer than sixteen 
points when any straight forward algorithm to compute min-circle is applied. This 
algorithm runs in linear time. 



As we will show later in the next section, the method of computation of intersection 
radius of a set of lines, , also uses the same principles. 



2.4 Intersection Radius Problem for Lines, Rays 
and Line Segments in the Plane 

In the previous section we presented a linear time algorithm to determine min-circle 
of a set of points P. In this section we shall design an algorithm to determine the 
minimum stabbing disk, min-disk , which intersects a given set of line segments, 
rays and lines. 

The real difficulty to apply the same technique to compute the intersection radius of 
line segments/ rays is that the bisectors of these can be curved linos. The nature of 
these depend on the relative positions of the line segments/rays in the plane. These 
may contain parts of bisector of two endpoints (a straight line), parts of bisector 
of an endpoint and a supporting line (a parabola) and a part of bisector of the 
supporting lines of both (an orthogonal pair of straight lines). It is the presence 
of parabolas in the set that does not allow us to apply the pruning technique to 
these bisectors. However, the curved parts, namely parabolas, in the bisectors of 
line segments owe their existence to the bisectors of line and points only. So, we 
can alternatively look at a similar problem in order to circumvent this difficulty. 
This problem is How can we determine the intersection radius of a mixed set 
of objects in the plane containing both lines and points. An added interest in this 
latter problem is that its solution is also used in the solution of the former problem. 

Further, each line segment is considered to be equivalent to two oppositely directed 
rays having the endpoints of the line segments as their respective tails. Thus the 
intersection radius problem for a set of /V/ 3 line segments can be reduced to the 
intersection radius problem for a set of 2/V/ a rays in the plane. The latter problem is 
solved by combining the prune and search strategy of Megiddo [21] with the novel 
idea of replacing a ray by a line or a point in the pruning step. This is done by 
breaking up the problem into the following two subproblems: 

Subproblem 4 Given a set of lines and points in the plane, compute the smallest 
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radius disk that intersects these. 

Subproblem 5 Given a set of rays, points and lines in the plane, show how a 
fraction of the rays can be replaced by lines or points such that the intersection 
radius of the new set is the sawe as that of the original one. 

We will use the subscripts /,;> and r to denote functions related to lines, points and 
rays, symbols L, P and R to denote a line, a point and a ray and denote finite sets of 
these by C,P and K, respectively. The functions f//(.r, ?y), <y ; ,(.r, ?/) and g r (x,y) have 
the following definitions: 



9i(x,y) = max<f/(L,(z,y)) | L 
g p (x,y) = max^ p (P,(j,T/)) P 



g r (x,y) = 



R 



where d/(L, (z, y)), d p (P, (a*, y)) and d r (R, (x, y)) respectively denote the (Euclidean) 
distance of the point (.r,y) from a line L, a point P ami a ray R. We define g(x,y) 



It can be easily shown that the functions gi(x,y),g p (x,y) and g r (x,y) are all convex 
so that ^f(x,y) is also convex. 

Let S be any finite set of lines, points and rays. The convexity of this last function 
enables us to answer two key questions about the minimum radius disk that inter- 
sects S. The first is that if we constrain the centre of the minimum radius stabbing 
disk to lie on a fixed line L then to which side of a given point (a,0) on L does the 
constrained centre lie. The second is that given a line L to which side of it does the 
centre of the minimum radius stabbing disk lie. Without any loss of generality, we 
can take L to be the x-axis of an orthogonal frame of reference. 



Let us answer the first question. Clearly, we can compute </(a,0) for the set S in 
O(N$) time, where N$ is the si/e of S. Let S' be the subset of objects whose 
distance from (a,0) is <?(a,0) or, more simply, those that touch the disk. Since 
<7(a,0) is convex, if the contact points of all the objects in S' lie to the left(right) of 
the vertical line through (rv,0) then the centre of the constrained minimum radius 
stabbing disk lies to the left(right) of (a,0). Otherwise, (a,0) itself is the required 
centre. 

The second question is also easily answered. We compute a minimum radius stabbing 
disk whose centre is constrained to lie on the line L. If the contact points of the 
objects in S' span an arc greatei than or equal to a semi-circle of the disk-boundary, 
then the computed disk is the required minimum radius spanning disk. Otherwise, 
the centre lies in the same halfspace, determined by L, as in which the mid-point of 
the chord of the above spanning arc lies. This again follows from the convexity of 
the function g(x,y). 

We first present a linear time algorithm to determine a minimum radius stabbing 
disk that intersects a given set of points and lines. When the objects to be intersected 
are only points, the algorithm of Dyer or Megiddo [11, 21] for the 1-centre problem 
can be used. The situation is more complex when lines are also included in the set of 
objects to be intersected. Our algorithm uses the basic prune and search technique 
of Dyer or Megiddo, referred to above. We describe the technique for lines only. We 
will then argue at the end of section that the addition of points does not change the 
underlying concept. 



2.4.1 Intersection Radius Problem for Points and Lines in 
the Plane 

Here we first present a linear time algorithm to determine min-disk of a set of lines, 
. Let the size of C be TV/. The min-disk of C can be computed by transforming 
the problem to a linear programming problem and then solving it using any of the 
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known methods [3]. The algorithm given here uses a different approach and applies 
the basic prune and search technique of Dyer and Megiddo [11, 21]. 

As in previous section, we progressively solve three similar problems. These problems 
are First, computation of the stabbing disk of minimum radius when the centres 
of stabbing disks are at a fixed point l\ computation of minimum radius stabbing 
disk when the centre of stabbing disks lie on a fixed line L and lastly, computation 
of the unconstrained min-disk. 

The first of these problems is solved in linear time by computing the maximum of 
the distances from each line in to point P. The required min-disk has this distance 
as its radius. 

In the solution of the other two, as in the case of points, the corresponding bisectors 
of pairs of lines play an important role These bisectors are either orthogonal pair 
of angular bisectors (for pair of non parallel lines) or parallel midway line (for pair 
of parallel lines). The lines are dropped in the computation of intersection radius 
by the application of the following rule: 



// the centre of nun-disk is localised in one of the regions, defined by the 
bisector (s) of two lines in the set L, then the line that is nearer to the 
localised region can be dropped. 



To prune away a constant fraction of lines in by applying this rule we need to 
localise the centre of the min-disk with respect to the bisector(s) of a constant 
fraction of pairs of lines. The localisation of centre with respect to a single line is 
done by computing the function g(x, y) as mentioned above. We discuss next how the 
localisation of the centre of min-disk, with respect to the bisectors (angular /parallel 
bisectors) of a constant fraction of pairs of lines in , in linear time, is done. 
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I A Constrained Version of the Problem 

We discuss the problem of computing the min-disk for a set of lines, , such that 
the centre of this disk is constrained to lie on a line //, which is assumed here to be 
the x-axis of an orthogonal frame of reference. We also assume that C has at most 
two lines parallel to L, since any other lying between these two does not change the 
intersection radius. Furthermore, out of these two lines we can prune away the one 
that is closer to L. So there is no loss of generality in assuming that there is at most 
one line parallel to L. To the remaining N[(> (N( I)) lines we apply the prune 
and search strategy in the following way. 

We first identify a group of roughly half the line-pairs such that one of their bisectors 
does not intersect an interval on L which contains the centre of the constrained min- 
disk. This can be done in linear time as follows. 

We compute the median x m of tho intersections of all the lines with L. Thru we 
determine the value of g(x m ,Q) a ^ this point and use this to locate the constrained 
centre on L with respect to (z m ,0). If the constrained centre, denoted by (x*,0), 
lies to the right (left) of this median point, the required half consists of those lines 
whose intersections lie to the left (right) of the median. We label the lines L,, with 
1 < i < [N//2J It will suffice to discuss the case in which the constrained centre 
lies to the right. 

Consider the sot of lino- pairs (7/ 2 i-i, ^2.), with \ < i < [N'i/4\. 

For each parallel line-pair it is clear that we can prune away the one that, is closer 
to x m . The pruning mechanism is non-trivial for non-parallel pairs, since we have 
to do a finer location of the constrained centre. Of these, we can prune away the 
line closer to x m if one of the angular bisectors is parallel to lino L. 

Each non-parallel pair (Z^.-i,/^,) nas an associated pair of angle bisectors. The 
intersection with L of one of these does not lie between the intersections of the lines 
themselves (Fig. 4). Let d t be this intersection point. 
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Figure 4: <f, does not lie between the intersections of lines L 2 , and L 2 i-i- 
d m lies to the left of x*. 




Fignrr 5- d m lies to the right of x* 

As before, we locate the constrained centre with respect to the median d m of those 
intersections among these that lie to the right of x m . Either x* c/ m , or we have 
the following cases: 



Case 1: z* < d m 

Consider a d t that lies to the right of d m . Since x* lies to the right of all 
the L.'s, the angle bisector that intersects L at d t , its associated pair, 
together with the lines that these bisect give rise to the configuration 
of Fig. 5. Since x* lies as shown we can prune away one of the lines of 
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the pair, (L 2l - 1,^21)- It can be easily shown that we can prune roughly 
one-eighth of the lines we started with. 



Case 2: x* > d n 



It is easy to see that in this case also approximately one-eighth of the 
lines we started with ran he thrown away (Fig. 4). | 



The above pruning takes 0(Ni) time. We repeat this process until no more lines can 
be pruned, when the optimal solution can be obtained by some brute force method. 
It can be shown that no more than eight lines are left at this stage. If the initial set 
of lines contained a line parallel to /,, the radius of the constrained min-disk is the 
maximum of the optimum value obtained and the distance from this parallel line to 
L. The total running time of this algorithm is easily seen to be O(Ni). 

US also contains points, each pruning step is carried out in two substeps. In the 
first substep we prune points, followed by lines in the second or the other way round. 
To prune points, we first pair them arbitrarily, and compute the median point of the 
intersections of L with the bisectors of these pairs. Then we determine the min-disk, 
Z), centred at the median point, for the set S. Next, to determine on which side 
of the centre of D on L the constrained centre lies, we examine how the points of 
tangencies of lines of <S, that touch the disk D, and the points of <S, that lie on its 
circumference, are distributed with rrspc-ct to the vortical line through the centre of 
D. When we prune lines next, again points of S are ignored similarly in the steps 
in which we determine x m and d m . 

In summary, we note that when pruning objects of one kind, the objects of the other 
kind become transparent whenever we need to find a point on L to serve as the centre 
of stabbing disks for all the objects in the set currently under consideration. 

If NI and N p be the number of lines and points in S respectively, it is easy to see that 
we prune away at least one-eighth of the total number of objects, viz NI + N p . We 
repeat this process until we cannot prune any more objects. The constrained centre 
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Figure 6 One of the bisectors of L, and L. (/?' in the figure) does not intersect the interior of 
LI 

can then be determined by a brute force algorithm. The whole process requires 
linear time. 



The Unconstrained Centre Problem 



Assume that C is the set of lines in S. We compute the inin-disk of in the 
unconstrained case. As before, we will indicate later how to handle the addition of 
points. 

We pair up the lines in C arbitrarily and compute the angle bisectors of each pair. 
In the degenerate case of a pair of lines being parallel, the angle bisectors reduce to 
a single line parallel to and equidistant from the lines that make up the pair. When 
a pair of lines have distinct angle bisectors, these make up an associated pair. 

The following observations are crucial. Given a pair of intersecting linos, if we can 
locate the region containing the centre in a quadrant defined by the angle bisectors 
of the pair then we can prune away one of the lines. For a pair of parallel lines we 
can do the same if we can locate this region in a halfplane determined by their "zero 
angle bisector". We indicate below how we do this for a fraction of such pairs. 
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Consider first the subset of vertical bisectors. We locate the unconstrained centre 
with respect to a median vertical bisector L\. Assume that the centre of min-disk 
lies in the left halfspace J\ of this median line. Clearly about half of the vertical 
bisectors do not intersect J\. We now compute the median slope of the non-vertical 
bisectors and pair up arbitrarily in this subset a bisector which has slope greater 
than the median slope with one which has a smaller slope. The bisectors in each of 
these pairs necessarily intersect as they have unequal slopes. Next we compute the 
median of all the unpaired bisectors whose slopes are equal to the median slope. We 
again locate the unconstrained centre with respect to this median bisector L 2 which 
has median slope. Assume that the centre lies in the halfplane J 2 below this line. 
Clearly about half of the median slope bisectors that lie above LI do not intersect 
J 2 . The centre now lies in J\ H ,7 2 . To be able to prune any lines, we need to refine 
the location of the centre still further as in the worst case, we may neither have any 
vertical bisectors nor any whose slopes are equal to the median slope. Therefore, we 
do the following with respect to the above pairs of intersecting bisectors. 

We first locate the unconstrained centre with respect to a vertical line through the 
median of the x-coordinates of their intersections. We can assume, without any loss 
of generality, that the unconstrained centre lies to the left, of this line. We now 
project, parallel to the median slope, the intersection points that lie to right of this 
line onto the y-axis. We then locate the unconstrained centre with respect to a line, 
parallel to the median slope, passing through the median of these projections. Again 
there would be no loss in generality if we assume that the unconstrained centre lies 
below this line so that now it lies in the lower left quadrant ////, determined by this 
line and the earlier one. 

This ensures that at least a fourth of the pairs of bisectors in this class have their 
intersections in the upper right quadrant U R. Consider one such pair. Since the 
bisector with slope smaller than the median slope does not intersect LL, at least 
one-eighth of the bisectors whose slopes are not equal to the median slope do not 
intersect LL. We note that this argument does not depend upon which quadrant the 
unconstrained centre lies in. Thus, we have a set of bisectors which do not intersect 
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Figure 7: The region J = LL n Ji H J 2 . 

the region J = LL D 7, n J 2 in which the centre of min-disk lies (Fig 7). In the 
worst case even this set of bisectors does not enable us to prune any lines because 
there may be no associated pairs of bisectors or bisectors due to a pair of parallel 
lines in this set. Further, in the worst case, it may happen that all of the associated 
bisectors of this bunch intersect J. 

We need to do one more final refinement of the location of the unconstrained centre 
with respect to these associated bisectors that intersect J. If we repeat the above 
steps with the above associated set of bisectors we get a region J', containing the 
centre and a fraction of these bisectors which do not intersect it. Each of these, 
together with its associated pair from the earlier set of bisectors we found that do 
not intersect J, contain in one of its quadrants the region J fl J 1 . Region J fl J' 
contains the unconstrained centre. Therefore one of the lines whose angle bisectors 
these are can be pruned. At the same time we can prune away one of the lines of 
each parallel pair whose "zero angle bisector" does not intersect J n J'. 

Calculation shows that at least [W//64J of the lines are pruned away. 

We repeat this process until no more lines can be discarded. It can be shown that 
there are no more than 64 lines in at this stage and we use some brute force 



26 



Xm L 



\ 



\ 



\ \ 

Figure 8: Ray R is replaced by the supporting line of R\ R' is replaced by its tail point. 

method to compute tlir unconstrained centre. 

When points are also included in the set, S, of objects to he intersected, we go 
about the pruning step in exactly the same way as in the constrained case. In each 
such step we throw away a well determined fraction of the points and the lines. 
Repeating this process we get a linear time algorithm for the unconstrained centre 
problem for S. 

2.4.2 Intersection radius for rays 

The constrained problem for a set 7 of N T rays in the plane can be reduced to the 
problem of pruning for a set of rays, lines and points as detailed below. 

For each ray R % G 7, consider the line L, normal to it and passing through its 
tail. We compute the median x m of the intersections of these normals with the 
constraint line L and locate to which side of this median on L the constrained 
centre lies. Suppose the constrained centre lies to the right of this median. Then 
for each normal which intersects L to the left of this median point we replace the 
corresponding ray by a line or point according to the following replacement rule : 



// the ray and the median point lie in the same halfspace of the two 
halfspaccs determined by the normal then the ray is replaced by the line 
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which contains the ray. Otherwise, the ray is replaced by its tail point 
(Fig. 8). 

It is important to note that neither of these replacements changes the radius of the 
relevant stabbing disks of the original set of rays. Thus, if N r be the number of rays 
in K, at least [N r /2\ rays are replaced by either a point or a line and our new set 
of objects consists therefore of lines, points and rays. Next, all the points and lines 
are considered for pruning as described in the 2.4.1. We thus discard a fraction of 
the rays from further consideration. These two substeps are iterated until no more 
objects can be discarded, when any brute force method can be applied to compute 
the min-disk centred on /,. It is easy to check that the algorithm runs in linear time. 

To solve the unconstrained version of the intersection radius problem for a set of 
rays, we need to replace a fraction of the rays by points or lines in linear time. This 
can be done as follows. 

As in the constrained case, we start with the normals through the tails of the rays. 
Proceeding identically as in the case of the unconstrained problem for lines, we 
determine a region J which contains the unconstrained centre and is not intersected 
by at least one-eighth of these normals. We do not need to iterate twice as in the 
case of lines. Only one iteration suffices for the replacement of a fraction of rays. 
The ray corresponding to each of these normals can therefore be replaced by a line 
or a point according to a similar replacement rule. 

// the ray and localised region J lie in the same halfspace of the two 
halfspaces determined by the normal then the ray is replaced by the line 
which contains the ray. Otherwise, the ray is replaced by its lad point. 

Thus in a single iteration we replace about one-eighth of the rays by lines or points. 
Next we use the two-step pruning process on the set of lines and points generated so 
far to throw away a fraction of them. Repeating these two substeps on the modified 
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set of rays, lines and points, and then applying any brute force method when no 
pruning takes place we get a linear time algorithm for intersection problem of line 
segments in the plane. 



2.5 Concluding Remarks 

In this chapter we have described optimal algorithms for computing the smallest 
radius disk which intersects a set of line segments in the plane using a novel approach. 
It would be worth investigating whether a similar approach can be used for other 
kinds of stabbing problems. 



Chapter 3 

An Optimal Algorithm for the 
Intersection Radius of a Set of 
Convex Polygons 



In the last chapter it was shown how the minimum stabbing disks and intersec- 
tion radii for finite sets of points, lines, rays and line segments can be computed. 
These were computed in linear time by combining the prune and search strategy 
of Megiddo [21] with the strategy of replacing line segments and rays by lines or 
points [3]. In this chapter, we enlarge the scope of this technique by showing that 
it can also be used to compute the intersection radius of a finite set of convex 
polygons in linear time. Moreover, it is immaterial if the set also contains other 
types of geometric objects such as points, lines, rays, line segments, half planes and 
wedges. In fact, we will show how to handle such a mixed set of objects in a unified 
way. 

The bisectors of line segments, as we have seen in the previous chapter, are made 
up of lines and (parabolic) curves. Here too the bisectors of convex polygonal disks 
are formed by not only parabolic curves but also regions of non-zero width. Such a 
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caso aiises when a pair ol convex polygonal disks have non-empty intersection and 
we have these regions of non zero finite area included in the bisector(s) of these. We 
present a method in this chapter in which this non linearity of bisectors does not pose 
any problem and we compute the minimum stabbing disk of a finite set of convex 
polygonal disks in linear time. But before doing this we first formalise the notions 
of replacement, substitution and localisation. These concepts were implied when we 
used these in their specific contexts in the previous chapter but nevertheless, were 
not discussed there. 



3.1 Introduction 

Let C be a finite collection of objects in the Euclidean plane. The stabbing problem 
consists of computing an object (the stabber) which intersects each member of C, 
Typically, the stabber could be a line, a disk, a line segment etc. and C could be a 
collection of points, lines, line segments, rays, circles, polygons or any mix of these. 
A survey of some recent results is available in the paper by Houle et al. [15]. 

The intersection radius of a finite collection of geometrical objects is the radius of 
the smallest closed ball that intersects all the objects in the collection. A classical 
problem in this area is that of finding the intersection radius of a finite set of 
points in the plane, which is also known as the 1-centre problem [30, 31, 6]. The 
corresponding disk is called the minimum stabbing disk. It was shown by Megiddo, 
and subsequently by Dyer [22, 12], how this can be solved in linear time. More 
recently, Welzl [3-1] has given a randomised algorithm for this problem which runs 
in expected linear time. However, until an earlier paper by Bhattacharya et al. [3] 
no attempt was made to extend this to a more complex collection of objects than 
points or to a collection containing different kinds of objects. 

New attempts have been made recently to find more complicated stabbers for 
the stabbing problem, or to find the best stabber which optimises some measure 
defined on the class of stabbers in question. Goodrich and Snoeyink [14] presented 



31 



an O(nlogn) algorithm to find a convex polygon whose boundary intersects each 
of n parallel line segments. Happaport and Meijer [24] showed that a perimeter 
minimising convex polygon that intersects each of n parallel line segments can be 
found in O(nlogrc) time. They have also extended their result to a set of isothetic 
line segments. Mukhopadhyay et al. [27] have shown that for a set of parallel line 
segments an area minimising convex polygonal disk can also be found in O(nlogn) 
time. Bhattacharya et al. [2] gave an O(n log n) algorithm for computing the shortest 
line segment that intersects a set of n line segments in the plane. 

Bhattacharya et al.[3] showed that when C is a collection of line segments the 
intersection radius can be found by combining the pi urn 1 and search strategy of 
Megiddo [21] with the strategy of replacing line segments with points or lines. This 
was discussed in the last chapter. In this chapter we enlarge the scope of this 
technique by showing that the intersection radius can also be found in linear time 
when C is a collection of convex polygons. Really, it is immaterial if C also contains 
other geometric objects like lines, points, rays etc. We show how it is possible to 
treat such a mixed collection of objects in a unified way. 

The organisation of this chapter is as follows. Section 3.2 contains the necessary 
geometric and algorithmic preliminaries. In Section 3.3 we describe the algorithm 
and analyse it in the following section. The last section contains conclusions and 
directions for further research. 



3.2 Preliminaries 

In the rest of this chapter we will adopt the following notation scheme: collections 
of objects will be denoted by letters in script style such as C,,T,. . ., objects by 
capital letters such as /,, /> /?, . . ., and points by small letters p,q, r, . . . Let the size 
of a collection C be denoted by NC- 

The algorithm that we shall describe in the following section is based on three 
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important techniques. These are: replacement of a complex geometrical object by 
a set of simpler objects, localisation of the centre of minimum stabbing disk in a 
more restricted region and thereafter filtration of some of the objects which do not 
play a role in determining the centre of the minimum stabbing disk. Together, the 
latter two make up the prune and search technique that was first formalised by 
Megiddo [21]. 

We discuss these strategies, replacement, localisation and filtration, in detail in the 
next three sections. 

Let C be a collection of ?i objects in the pla.no where each object, is either a point, a 
line, a line segment, a ray, a wedge or a convex polygonal disk. The problem is to 
determine the minimum stabbing disk, min-disk, of the collection C. 

In the following discussion J is a fixed subset of plane containing the centre of min- 
disk. In the general case, ./ is initially the whole plane, and gets smaller as the 
computation proceeds. 

In the constrained case, J is always a subset of the line on which the centres of 
stabbing disks are constrained to lie. All the discussion in this section is also 
applicable to the constrained case where J is such restricted to a line. 

Let C and C denote an object and a collection of objects respectively. The distance 
of a point p from C, d(p, C), is the shortest Euclidean distance of p from C '. The 
distance of p from C, d(p,C), however, is the largest of all distances from p to the 
objects in C. The stabbing radius, $R(J,C), where the centre of the stabbing disks 
are constrained to lie in J , is the minimum of the distances from points in ,7 to C. 
The intersection radius, IR(C), is the unconstrained stabbing radius. In summary, 
we have 



d(p,C) = inf%, 9 ), 

(jfcC- 

d(p,C) = max d(p, C) 
SR(,7,C) = m\nd(p,C), 
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Figure 9- A polygonal disk, and its substitution by wedges (one of the wedges seen slightly shifted) 

IR(C) = SR( 2 ,C), 
where d(p, q] is the distance from p to q. 

The distance bisector of two objects S and T is defined as the set of all points for 
which the distance to 5 is equal to the distance to T, i.e., it is the set 



3.2.1 Object Replacement 

To simplify the computation it is helpful to divide a convex polygonal disk into a 
set of elementary parts such that the distance function is not modified. This allows 
us treat these in a uniform manner in the algorithm. To see how we can divide a 
convex polygonal disk, let us first see how the distance from a point p in 9R 2 to a 
Convex polygonal disk is computed. This distance is either equal to the distance 
from p to some vertex, or it is equal to the perpendicular distance to some side, 
or it is zero when the point lies inside the disk. A wedge is the non-reflex region 
bounded by its two infinite sides. It can be easily seen therefore that the distance 
from p to this disk is the maximum of the distances to the wedges formed by taking 
the vertices and their adjacent sides (Fig 9). So if a convex polygonal disk in C is 
substituted by these wedges then the intersection radius of the resulting collection 
does not change. 
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V 

Figure 10: A line segment, and its substitution by two rays (seen shifted) 

Each line segment can be similarly replaced by two oppositely directed rays whose 
intersection is the line segment itself (Fig 10). The distance of the line segment from 
a point p in & 2 is either equal to the distance of the endpoints or it is equal to the 
perpendicular distance to the supporting line; this distance is clearly equal to the 
maximum of the distances of p from the two rays. 

These examples motivate the following theorem. 

Theorem 3.1 Let C be an object and T> a collection of objects such that d(p,C) = 
d(p, T>) for all p J. Then for any collection C, 



Proof: From the definitions of intersection radius and distance functions, we have 



SR(J,CU{C}) = min 

= min m&x{d(p,C),d(p,T>)} 

= min d(p,C U T>) 

= SR(J,CUX>). 

This proves the theorem. | 

When C is a line segment, a ray, a convex polygonal disk, a halfplane or a wedge, 
we can apply Theorem 3.1 to leplwr C by some simpler objrct(s). 
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Ray is replaced . Ray it replaced 

by tail point by supporting line 

\ 

Figure 11: A ray R, its associated normal and regions 

In the two examples mentioned above J is the whole plane. This allows unconditional 
replacement of line segments and convex polygonal disks by rays and wedges, respec- 
tively. Let us take another set of examples, when J is only a proper subset of the 
plane, C is a ray, a halfplane or a wedge, and the replacement is done conditionally. 
Let R be a ray in C. The normal to this ray through its tail point divides the plane 
in two halfplanes. If we localise the centre of min-disk in one of these halfplanes 
then the distance from any point of the localised region to the ray is equal to either 
the distance to the tail point or the perpendicular distance to the supporting line 
of the ray R (Fig. 11). This means that we can replace R either by its tail point or 
its supporting line. 

Similarly if J is localised in the interior of a halfplane, //, in C then the distance 
from any point of J to H is zero, otherwise if J is localised in the interior of H then 
the distance from any point of J to H is equal to the distance to the boundary of 
H . We can respectively discard H or replace H by its boundary line in C if these 
cases arise (Fig 12). 

Let us draw an outward normal to each of the sides at the apex of a wedge W '. 
The whole plane is then divided into four (unequal) quadrants by these normals and 
sides of W. If we somehow localise the region J in one of these quadrants then the 
distance of the points in J from the wedge has one of the following values: 0, when 
J is localised inside the wedge; perpendicular distance to one of the sides, when 
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Replacement by 
Boundary Line 



Deletion of i_ 1 

Half Plane 



Figure 12: A halfplane //, its associated boundary line and regions 




Deletion of 
Wedge 

y 

Wedge 

Figure 13: A wedge W , its associated normals and regions 

J is localised in the region bounded by a side and the normal to it; or, distance 
to the apex, when J is localised in the region bounded by two normals. We can 
respectively discard W , replace W by the relevant side, or replace W by the apex, 
in these cases(Fig 13). 

We shall henceforth view a convex m-gon as the collection of m wedges defined by 
the vertices and the sides incident on them. Our problem, therefore, is equivalent to 
that of computing a min-disk for a set of A r w(= m > sum taken over all the convex 
polygons) wedges. Likewise, we shall view a line segment as a pair of oppositely 
directed rays defined by its end points [3]. 
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Figure 14- Stabbing disk is a rain-disk 

3.2.2 Localisation of the Centre of M in- Disk 

In the discussion below we assume that the objects are convex sets of points so that 
all distance functions are convex. 

An object C is said to touch a disk of radius r centred at point p if r = <f(p, C) , 
and the vectors from p to the points of contact are said to be its contact vectors. 
The following theorem characterises min-disk. 

Theorem 3.2 A stabbing disk of radius r centred at point p is min-disk iff all of 
its k contact vectors, r,, with 1 < t < k, are linearly dependent satisfying 



where A, > 0, for 1 < t i < k, with some A, ^ and k > 2. 

Proof: The proof makes use of the behavior of the distance function in the 
neighbourhood of p, which depends only on the contact vectors of the stabbing disk. 

Let the contact vectors r, satisfy the relationship given in the theorem. Then for 
any arbitrary vector v, we have 
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and hence, (r, v) < for some t, with 1 < t < k. Let this object be C. Now, 
if p is displaced to p' by a small amount 6v in the direction of v and if 6v' is the 
corresponding displacement of the new contact point, p' t , of the object C (Fig 14) 
then the square of the new distance, <f(, of the new centre p' from the object C is 



= |r.| s + \6v - *vf - 2(r, <v) + 2(r, <5v') 



Since C is convex and r, is normal to it, r, 6v' is non-negative and this justifies 
the inequality above. Since d(p,C) is the maximum of all distances from p to the 
objects in the collection C, d(.,C) will also increase in the direction v. Thus 



> r 

> d(p,C). 

Hence, as v is an arbitrary vector, p is a local minimum of the function d(.,C). Since 
<f(.,C) is convex, p is a global minimum too. 

For the proof of the converse, we assume that the said disk is the min-disk. Then we 
prove that the centre p lies in the relative interior of the convex hull of the contact 
points. 

We will prove this by contradiction. We assume the contrary that the centre of min- 
disk lies does not lie inside the convex hull of contact points. Then we can compute 
a line which separates the point p and the convex hull. Let v be the normal vector 
to this line, contained in the same halfplane as the convex hull. Then v r t > for 
all i, with I <i < k. From this it can be seen that d(p,C) strictly decreases along v, 
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Figure 15: Centre of Min-Disk lies inside convex hull of contact points 

< Ir.-H 2 

< |r,| 2 - fy (2r, - v) 

< |r,f 

< r 



This means that the radius of the min-disk is not minimum. Since this cannot be, 
therefore the assumption that the centre of the min-disk lies outside the convex hull 
of contact points is incorrect. 

Let the vectors to p and the contact points be p and p t , with 1 < i < k, respectively. 
Now, an interior point of a convex hull can be written as a positive linear combination 
of the extreme points. Thus 
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where A, > 0, with 1 < t < k. Hence, 



A,( Pt -p) 

t<k 



and therefore 



The above theorem is also applicable in the constrained case with a slight modifi- 
cation that here the contact vextors projected onto the constraint line are positive 
linearly dependent. Once we have computed a stabbing disk, if it is not a minimum 
one as determined by Theorem 3.2, we can further localise the region in which the 
centre of min-disk lies with the help of the theorem below. 



Theorem 3.3 Given a stabbing disk of C centred at p, with k contact vectors r l; 
with 1 < i i < k, the centre of min-disk lies in the set J' given by 



i{n^}. 



where H t is the half space, normal to r,, that passes through p and contains the 
contact point and the tangential object. Furthermore, 



Proof: Let the radius of stabbing disk be r. We have to show that centre of 
min-disk lies in every H t . Let the corresponding object to Hi be C. The radius of 
min-disk has to be smaller than or equal to r. Therefore the centre of min-disk is 
at a distance smaller than or equal to r from every object in C and in particular 
from C. The feasible region of this centre will thus be a subset of //, (because C 
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Figure 16: Localisation of Centre of Min-Disk 

is convex). Every point not in //, is at a distance larger than r from the object C 
(Fig 16). This proves the theorem. | 

We can use this theorem to localise the centre of min-disk with respect to any 
arbitrary fixed line, Z,, in the plane. We first compute the minimum stabbing disk 
such that its centre is constrained to lie on L and then compute the region given 
by the expression of the theorem with this disk as the reference. Since the region 
contains only those points where the value of stabbing radius is smaller than the 
radius of given disk, therefore it necessarily does not contain any point of L. Thus 
we get a region that is fully contained in one of the halfplanes determined by L. 
Further, if we localise the centre of min-disk to a region J' we then onwards replace 
and filter objects with respect to newly located region J' and not the earlier located 
region J. 



3.2.3 Pruning or Filtering Objects in C 

We now need a suitable criterion to discard objects from the collection C that are 
irrelevant in determining the centre of min-disk. The following theorem provides 
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Figure 17: Filtering lines and points 



such a criterion. 



Theorem 3.4 Let there be a pair of objects in C such that region J is contained in 
one of the regions determined by their distance bisector ($). Then one of these objects, 
whichever is nearer to J, can be discarded from C without affecting the intersection 
radius. 



Proof: From the definition of the distance bisector, if the region in J is contained 
in one of the regions defined by the bisector then, for every point p in 7, distance of 
p from one of the objects is always smaller than the distance from the other object. 
Let this nearer object be C. Since the intersection radius is the maximum of the 
distances of objects of C over the points of the region .7, the object C , its distance 
being smaller than the distance of the other object in the pair, does not play a 
role in the determination of min-disk. Hence, the object C can be deleted from the 
collection C. | 



For example, refer to Fig. 17 above. The angular bisectors of two lines LI and LI 
and the perpendicular bisector of two points PI and PI &r e shown. The nearer of the 
lines, LI, is deleted because region J (the shaded region in the figure) is contained 
in the one of the four quadrants defined by the angular bisectors of LI and L^- 
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Similarly, the nearer of the points, P^ is deleted because the region J is contained 
in one of the halfplanes defined by the bisector of PI and P 2 . 



3.3 Intersection Radius Problem for Convex Poly- 
gons 



In this section we present a linear time algorithm for computing the min-disk which 
intersects a collection C of convex polygons. In fact, we will be more general and 
assume that C also contains points, lines, rays, line segments, halfplanes and wedges. 
Since each line segment can be considered to be equivalent to two oppositely directed 
rays having the end points of the line segments as their respective tails, and a convex 
polygon as a collection of as many wedges as vertices (refer 3.2.1), we will assume 
that C consists of lines, points, rays, halfplanes and wedges only. 

We first solve a constrained version of the problem and then the unconstrained one. 
The motivation behind this is as follows. The intersection radius is the minimum 
of all stabbing radii over the domain of the plane. This problem can be formulated 
as a non- linear programming problem. The usual method of solution of such a 
problem in E d consists of locating the optimal point with respect to a hyperplane 
in one dimension lower, by first solving the problem restricted to this hyperplane, 
and then locating the optimal point by computing the gradient of the minimising 
function. 

Furthermore, to replace or filter objects, we need to locate the centre of min-disk 
with respect to either bisectors of pairs of points, bisectors of pairs of lines, the 
wedges or the normals to rays. So we apply the same technique, as in non-linear 
programming, to localise the region containing the centre. We note that we need 
not do this for every object in C in a single step. It suffices to do this only for a 
fraction of these objects. So, we successively localise the centre of min-disk with 
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Figure 18: Localisation of Centre of Min-Disk in Linear Time 

respect to a group of lines, which are chosen in such a way that there is always a 
fraction of the objects in C which are either discarded or replaced. 

First we provide an efficient linear time implementation of the construction suggested 
by Theorem 3.3. Then we solve the constrained and unconstrained cases of the 
intersection radius problem in successive subsections. 



3.3.1 An Efficient Implementation of Theorem 3.3 



Suppose we are given a stabbing disk centred at c. We have to locate the centre of 
min-disk in a region, J, such that c lies on its boundary. We do this in the following 
way. 

If the stabbing disk does not touch any object then its radius can be shrunk until 
it touches at least one object. This can be done in linear time by computing 
the maximum of the distance from c to the objects in C. Suppose there are k 
contact vectors of the stabbing disk with k > 1. If we compute the intersection of k 
halfplanes as in Theorem 3.3 by constructing the convex hull in the dual plane, then 
we will need 0(k log k) time. This can take 0(Nc log NC) time in the worst case. In 
our application, however, this intersection can be computed in linear time because 
all these halfplanes have point c on their respective boundaries. Computation of 
their intersection is same as the convex hull computation in the dual plane of a 
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set of collinear points. We only need to compute the two extreme halfplanes and 
the optimal centre lies in the intersection of these two halfplanes which is a wedge 
(Fig 18). If the intersection of the halfplanes is a single point, namely c, then c 
is the centre of min-disk. In the constrained case we need to locate the centre of 
minimum stabbing disk on a line L. For this we compute the part of the line L that 
lies inside this wedge. 



3.3.2 The Constrained Centre Problem 

In the constrained version of the problem we determine a min-disk whose centre 
is constrained to lie on a line L. Let P, C, 71, Ti. and W be disjoint subsets of 
C containing points, lines, rays, halfplanes and wedges respectively such that their 
union is C. We show how we can process points, lines, rays, halfplanes and wedges 
separately. We then unify these into a single iteration of the algorithm. 

Points: We pair up the points in P arbitrarily and compute the median of the 
intersections of their bisectors with the line L. The centre is then localised 
with respect to this point on one of the half lines. We filter out at least one 
fourth of the points corresponding to the bisectors that do not intersect the 
localised region. We repeat this step until no filtering takes place. 

Lines: We pair up the lines in C arbitrarily and compute their distance bisectors 
(angular bisectors for non-parallel lines, and mid- way line for parallel lines). 
We divide the intersections of these bisectors with the line L equally into four 
intervals on L. The centre of min-disk is localised in one of the intervals by 
doing a binary search on the boundary points. We filter at least one fourth 
of the lines corresponding to the distance bisectors that do not intersect the 
localised region. We repeat this step until no filtering takes place. 

With rays, halfplanes and wedges we proceed in a slightly different way. 
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\ 

Figure 19: Intersection of objects with J in the constrained problem 

Rays: We draw normals at the tail points of the rays in 7 and compute the median 
of the intersections of these normals with the line L. The centre of min-disk is 
then localised with respect to this point. We can replace a ray, whose normal 
does not intersect the localised region, either by its tail point or its supporting 
line. We replace at least one half of the rays in each iteration. We repeat this 
step until no replacement takes place. 

Halfplanes: We localise the centre of min-disk on the line L with respect to the 
median of the intersection of line L with the boundaries of halfplanes in H. 
We replace a halfplane in *H by its boundary line or discard it if its boundary 
line does not intersect the localised region. We replace/filter at least one half 
of the halfplanes in each iteration. We iterate until no filtration/replacement 
is done. 

Wedges: For each wedge in W, we draw outward normals to the sides at its apex. 
We divide the intersection(s) of the normals and the sides with the line L 
equally into four intervals of L. We then localise the centre of min-disk to 
one of these intervals. All those wedges (at least one fourth) whose sides and 
normals do not intersect the localised region are either filtered out or replaced 
by a line or point. We iterate until there is no replacement or deletion. 

It is easy to design an O(NclogNc) algorithm, based on the above facts. We 
first convert all the rays in C to lines or points, then we filter or replace all the 
halfplanes and wedges in C, and finally compute the intersection radius of a collection 
containing only lines and points. To obtain a linear time algorithm, we have to treat 
the objects in a more unified way in each iteration of the algorithm. 
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The basic idea is as follows: The derived lines obtained as above from different 
objects of C are classified as either active or inactive, depending on the fact, whether 
region J which contains the centre of min-disk is localised with respect to these. 
Further, active lines among these are assigned following weights: derived lines of 
wedges with four and three intersections with J - 16; derived lines of wedges and 
line pairs with two intersections with J - 18; derived lines of rays and line- pairs with 
one intersection with J - 27; and derived lines of point-pairs with one intersection 
with J - 36. At the beginning of computation all the derived lines that intersect 
L are active. In each iteration we first localise the centre of min-disk on L with 
respect to the weighted median of intersections of line L with the active derived 
lines. The locating of the centre of min-disk is done by invoking Theorem 3.2 and 
Theorem 3.3 to determine on which side of the weighted median intersection point 
the constrained centre lies. Finally, we replace and filter those objects, none of whose 
derived lines are active, by invoking Theorem 3.1 and Theorem 3.4. We repeat this 
until no active derived lines are made inactive. 

It can be easily seen that in every iteration half of the active lines of about half 
the total weight are inactivated. New derived lines of at most only three-fourth of 
this discarded weight are added. So there in net reduction of weight by a constant 
fraction. Hence the algorithm is linear in total weight. 

The details are given in the algorithm below. T> is the active derived set of lines. 

Algorithm 1 CONSTRAINED-CENTRE 

Input: Bisectors of point-pairs and line-pairs, normals to rays, 
boundary lines of half planes and sides of and normals to wedges. 
Output: Centre of constrained min-disk 
begin 

V < input set 
do 

Compute the weighted median intersection, P, of T> with L 
Locate J with respect to P on L that contains the 
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constrained centre 
Update V 

Filter and replace relevant objects of the collection C 
Update 

while X> is updated 
enddo 

Determine the constrained centre by some brute force method 
end | 

The proof of correctness of the algorithm CONSTRAINED-CENTRE is a direct conse- 
quence of the theorems in Section 3.2. 



3.3.3 The Unconstrained Centre Problem 

The tasks of filtration and replacement are more involved in the unconstrained case. 
We first describe a method by which we can localise the region that contains the 
centre of min-disk which does not intersects a finite fraction of a given set of lines. 
As -in the constrained case, we then show how we can process points, lines, rays, 
halfplanes and wedges separately. Finally, we unify these separate steps into a single 
iteration of the algorithm. 

Let LOCALISE be the procedure which, given a set of lines L with integral weights, 
determines a region J that contains the centre and is disjoint from at least one 
eighth of the total weight of lines. Also, our algorithm terminates if the centre of 
min-disk is found during a call to the procedure CONSTRAINED-CENTRE in LOCALISE. 
So we only need to consider the case in which the algorithm does not terminate in 
this manner. 

We first compute the weighted median slope, s m , of the non-vertical lines in and 
then divide the lines into three sets, >, < and -, having slope greater, smaller 
than and equal to s m respectively. Let be the set of vertical lines in . We pair 
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the lines in > and < arbitrarily using their weights. We then compute intersection 
points of these pairs. The weights of both members of the pair should be equal and 
this weight is assigned to the intersection point. For example, a line with weight 3 
can be paired twice once with a line of weight 2 and second time with a line of 
weight 1 whereas the intersection points have weights 2 and 1 respectively. We can 
ensure the number of pairs to be linear as all the lines have integral weights only. 
Note that if the total weights of the sets > and < are unequal then some of the 
lines may remain unpaired. 

Next we compute a vertical line, L y , which divides the total weight of the intersection 
points of the above pairs and the vertical lines in C v into two equal halves. We 
localise the centre of min-disk in one of the halfplanes determined by Z/ v , by using 
the solution of the constrained centre problem with L y as the constraint line. Let 
J y be the halfplane determined by L y in which the centre of min-disk lies. 

We then compute a line L m with slope s m which divides the total weight of points 
and lines intersections points of the above pairs which lie outside J y , and lines of 
= into two equal halves. We again locate the constrained centre on this line and 
determine on which side of it the centre of min-disk lies. Let J m be this halfplane 
so that the centre of min-disk now lies in J v fl J m . 

It is easy to see that at least one eighth of the total weight of lines in C are disjoint 
from J y fl J m . These lines cross the opposite quadrant J v C\J m . 

The description of LOCALISE is as follows: 

Algorithm 2 LOCALISE 

Input: Set of Lines C 
Output: Localised Region J 
begin 

Compute the intersections of an arbitrarily (equal weighing) 
paired set of lines in , one of larger slope 
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than the weighted median slope and the other of smaller slope 

Localise Centre in J v with respect to a vertical line which 
halves the total weight of intersection points and 
vertical lines in 

Localise Centre in J m with respect to a line of weighted median 
slope which further halves the total weight of intersection 
points in 7~ v and weighted median slope lines in C. 

J +J m r\J y 

Return J 
end | 

Once we have localised the centre of min-disk, we repeatedly apply Theorem 3.1 
and Theorem 3.4 on wedges, halfplanes, rays, lines and points to prune/replace the 
objects in the collection C. We discuss separately the cases of points, lines, rays, 
halfplanes and wedges. 

Let P, , 7, H and W be the subsets of C containing all points, lines, rays and 
wedges in C respectively. 

Points: We pair up the points in P arbitrarily and compute the perpendicular 
bisectors of these pairs. We use procedure LOCALISE with this set of bisectors 
and filter one point corresponding to every bisector that does not intersect 
the localised region. Per iteration this filters out at least l /\e of the points 
(Fig. 20(a): The centre is localised in Lower Left quadrant, LL and the point- 
bisector crosses Upper Right quadrant, UK). 

Lines: Lines are filtered by localising the region containing the centre of min-disk 
in a quadrant defined by the bisector(s) of a pair of lines. When the lines are 
parallel the distance bisector is a line parallel to and equidistant from them. 
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We will call a pair of angular bisectors mates of each other in the following 
discussion. When we invoke procedure LOCALISE with this set of bisectors as 
input, there is no way to ensure that there exists at least a fraction of these 
bisector pairs which do not intersect LL. So, we have to do the localisation 
twice in the following way. 

We first pair the set of lines in arbitrarily and compute the bisector(s) 
of these pairs. We localise the region containing the centre of min-disk with 
respect to this set of bisectors. Let LL be this localised region. We consider all 
the bisectors that do not intersect LL. In the worst case, all these are angular 
bisectors such that their mates intersect LL. We invoke LOCALISE again with 
these intersecting mates as input. A fraction of these mate-bisectors do not 
intersect the localised region, say LL', returned by this second invocation of 
LOCALISE. These together with there mates from the previous invocation of 
LOCALISE do not intersect region LL fl LL' which contains the centre of min- 
disk. We filter one of the lines of each pair of lines whose distance bisector(s) 
does not intersect LL D LL' . This filters out V<54 of the lines in C (Fig. 20(b)). 

Rays: Rays in K are replaced by localising the region containing the centre of min- 
disk in one of the halfplanes defined by the normal at the tail point of a ray. 
We do this by invoking LOCALISE with the set of normals at tail points of rays 
as input. This replaces at least i/ 8 of the rays in ft (Fig. 20(c)). 

Halfplanes: Halfplanes in 7i are filtered/replaced by their boundary lines by lo- 
calising the region containing the centre of min-disk either in the interior of 
these or in the interior of their complement. We do this by invoking LOCALISE 
with the set of boundary lines. In each iteration at least l /s of the halfplanes 
are filtered/replaced in this way. 

Wedges: Wedges in W are filtered /replaced by localising the region containing the 
centre of min-disk in one of the four quadrants at the apex of a wedge (Fig. 20(d)). 
We then replace or filter the wedge, as the case might be. We call the four 
lines, two sides and two normals, associated with each wedge to be mates of 
one another. 
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Figure 20: Localisation of unconstrained centre in LL with respect to derived lines 



53 



We proceed exactly as in the case of line bisectors by first invoking LOCALISE 
with the set of derived lines of the wedges as input. Since for crossing derived 
line in the UR quadrant at most three derived lines of the corresponding 
wedge can intersect LL, we invoke LOCALISE twice. The first call is with the 
line/middle of these intersecting lines. Let LL' be the region returned. We 
invoke LOCALISE again with the remaining line, if any, that intersects LL' and 
whose mates do not. Let LL" be the region that is output. Thus we have a 
fraction of wedges in W whose associated lines do not intersect LLC\LL'r\LL". 
This replaces/filters at least V 2 56 of the wedges. 

We can get an O(Nc log NC) algorithm for computing the min-disk if we first convert 
all the wedges, halfplanes and rays to lines/points as above and then compute the 
min-disk of the resulting set consisting of lines and points only. To obtain a linear 
time algorithm, we do not differentiate between the bisectors of points and lines, 
associated lines of rays, halfplanes and wedges in the invocations of LOCALISE. 

We classify the derived lines as active or inactive in a similar way as in the con- 
strained case. We assign weights to the active derived line in exactly same way as 
the constrained case. At the beginning all the derived lines are active. We invoke 
LOCALISE with the set of active derived lines and if any active derived line does not 
intersect the localised region then we make it inactive for the rest of the computation. 
We replace/filter those objects in C none of whose derived lines are active. 

We repeat these steps until no update occurs. We later show that this algorithm 
runs in linear time. 

We give a description of the algorithm for the unconstrained centre problem below. 
We will denote the active derived lines by >. 

Algorithm 3 UNCONSTRAINED-CENTRE 

Input: Bisectors of point-pairs and line-pairs, normals to rays, 
boundary lines of halfplanes, and normals to and sides of wedges 
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Output: Min-disk 
begin 

V < input set 
do 

LOCALISE with V 

Update D 

Filter/Replace objects in C 
while there is some update 
enddo 

Determine the centre by some brute force method 
end | 



3.4 Analysis of the Algorithm 



In this section we establish that the algorithm runs in linear time. First we will 
analyse the constrained case and then the unconstrained one. 

In every iteration of the algorithm some of the wedges, halfplanes, points and lines 
are dropped or some of the wedges, halfplanes and rays are replaced by points or 
lines. The DAG above (Fig. 22) illustrates this conversion. Further each localisation 
of J with respect to some derived line reduces the weight by at least one fourth. 
For example, when a wedge having four intersections with J is converted to another 
wedge with three intersections with J then the net reduction in weight is 16 
from 64 to 48. We show this reduction in Fig. 21. The number of intersections are 
denoted in the subscripts, and the number of objects resulting from the conversion 
is written on top of arrows. Line pairs and point pairs are denoted by CP and PP 
respectively. 

Let NU denote the number of objects of type U. 
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Figure 21: Reduction of Weights Per Iteration 





Figure 22: Conversion among Object Types 
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Theorem 3.5 Algorithm CONSTRAINED-CENTRE is linear in NT + N c + NK + N* + 

Proof: 

Let WT> be the total weight of active derived lines. 

At least one of the lines derived from an object or a combination of objects in 
C intersects L. Otherwise we can straightaway filter or replace the corresponding 
objects. Then we have following bounds on Wp 



+ 277V /2 + 
+ 277V* 



Therefore, 

wi>/48< NP + NC + NK + NH + NW < 2^/27. 

Hence 

WT> = e(Np + NC + NK + NK + N w ). 

In every iteration the derived active lines of at least half the total weight are made 
inactive. The new derived lines of at most three fourth of this weight are added. 
So at least a fraction (around one eighth) of total weight is reduced. The pruning 
step is also linear in the number of objects which has the same complexity as the 
number of active derived lines. Therefore the algorithm is linear in Wp. Hence the 
algorithm CONSTRAINED-CENTRE is linear in Np + NC + NK + NH + Ww as claimed. 



The procedure LOCALISE, which locates the unconstrained centre with respect to 
constant fraction of weighted lines, uses the algorithm CONSTRAINED-CENTRE as its 
basic routine. So, linearity of latter implies the linearity of former. We now prove a 
linear time bound for the algorithm UNCONSTRAINED-CENTRE. 



57 



Theorem 3.6 Algorithm UNCONSTRAINED-CENTRE is linear in JV W + N H + N* + 



Proof: 

Let WT> be the number of active derived lines in the algorithm UNCONSTRAINED-CENTRE. 
Similar to the constrained case above, at least one of the derived lines of each object 
or the combination of objects in C is active. Otherwise we filter or replace the 
corresponding object/objects. Thus we have following bounds on Wp 



Therefore, 

W V /M< NP + NC + NK + NH + NW < 2^/27, 

and hence, 

W v = e(N P + N c + N* + TV* -f N w ). 

The pruning step is linear in the number of objects in C. Furthermore, in every 
iteration, at least one eighth of the total weight of derived active lines is made 
inactive and new active lines of at most 3 / 4 weights are added, with the consequence 
that the total weight is reduced by 1/32- Therefore the algorithm is linear in HV 

Hence the algorithm CONSTRAINED-CENTRE is linear in TVyy + NH 4- NK + N c + Np 
as claimed. | 



3.5 Concluding Remarks 

In this chapter, we have described an optimal algorithm for computing the smallest 
disk that intersects a finite collection of geometrical objects, containing convex 
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polygonal disks, wedges, halfplanes, lines, points, line segments and rays by com- 
bining the prune and search technique of Megiddo with the novel idea of replacing 
complex geometrical objects by simpler ones. 

It would be worth investigating whether this approach can be used to compute the 
smallest intersecting disk of a collection of objects that includes simple polygons 
and also look at higher dimensional generalisations of the problems studied here. 

In the next chapter we again see one more application of the modified prune and 
search, which was used in this chapter, in the computation of a centre point for a 
finite planar set of points. 



Chapter 4 

Computing a Centrepoint of a 
Finite Planar Set of Points in 
Linear Time 



The notion of a centrepoint of a finite set of points in two and higher dimensions is 
a generalisation of the concept of the median of a set of reals. In this chapter, we 
present a linear time algorithm for computing a centrepoint of a set of n points in the 
plane, which is optimal compared to the 0(nlog 3 n) complexity of the previously 
best known algorithm. We use suitable modifications of the ham-sandwich cut 
algorithm in [23] and the prune and search technique of Megiddo [21] to achieve this 
improvement. 



4.1 Introduction 

We all have an intuitive idea as to what phrases like "the very centre of the square" 
or "the very centre of the city" mean. To capture this intuition in a quantitative 
way, the centre of a set of n points, P, in ^ is defined as the maximal subset of 3^ 
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such that any closed halfspace intersecting this subset contains at least \n/(d + 1)] 
points of P [35]. This subset is non-empty for any finite configuration of points (see, 
for example, [13]). Furthermore, it is closed and convex. A centrepoint is a member 
of the centre of P. 

On the real line 3fc, a centrepoint is no other than a median of P. Thus a centrepoint 
can be viewed as a generalisation of the median of a set of reals. On the other hand, 
the centre can also be viewed as a particular k-hull of P. The k-hull oi Pisa maximal 
subset (closed and convex) of $?* such that any closed half space intersecting this 
subset contains at least k points of P. For instance, the 1-hull of P is its convex hull 
and the centre is its \n/(d + l)~|-hull. The property of balanced partitioning makes 
the centre point useful for efficient divide and conquer algorithms in geometrical 
computing and large scale scientific computing[26, 32, 25, 36]. Recently Donoho 
and Gasko have suggested that centre point can be used as "robust" and high 
"breakdown point" estimators for multivariate datasets [10]. 

The interesting algorithmic problem of computing a centrepoint has been considered 
by various researchers. Cole et al gave an 0(nlog 5 n) algorithm for computing a 
centrepoint of a planar set of points (9). Subsequently, Cole improved this bound 
to 0(nlog 3 n), using the powerful technique of slowing down a sorting network 
[8]. In this chapter, we propose an optimal linear time algorithm for computing a 
centrepoint of a planar set of points by using suitable modifications of the ham- 
sandwich cut algorithm for a pair of separable point sets [23] and the prune and 
search technique of Megiddo [21 J. 

Linear time algorithms, however, were known for computing an approximate or e- 
centrepoint [19, 32, 23]. We obtain this weaker type of centrepoint if we decrease 
the lower bound, in the above definition of the centre, to fn(l - e)/(rf+ 1)1, where 
< e < 1. Actually, Megiddo [23] only gave an algorithm for computing a partition 
of a (planar) set of n points with two lines such that each closed quadrant contains at 
least [n/4j points. An algorithm for computing an e-centrepoint, where < e < 1/4, 
is implicit in this. 
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The thesis proposes an optimal algorithm for computing a centrepoint of a planar 
set of points by using an interesting modification of Megiddo's prune and search 
technique[21]. This consists of adding a few extra points in each pruning step so that 
a subspace of the original solution space is retained, while ensuring a net deletion 
of points. In the description of our algorithm, we assume the usual RAM model of 
computation; the point set P, however, is not assumed to be in general position. 

This chapter is organised as follows. In section 4.2, we discuss which points to prune. 
In section 4.3, we describe the method used to find these points. The algorithm is 
presented in section 4.4. Section 4.5 contains an analysis of the time complexity of 
the algorithm. Concluding remarks are given in section 4.6. 



4.2 What to Prune 

Let P be a finite set of points in the plane. In the subsequent discussion we use 
the following notations. We denote the centre of P by CEITRE(P) and the fc-hull of 
P by WLL(k,P). We use the notations P H , PGH, PFGH, to denote the points 
of P contained in //, G D H, F D G fl #, . . . respectively, where F, <7, H, . . . are 
any closed or open halfplanes. We denote the complement of a set 5 by 5. As we 
frequently need to use the numbers \\P\/$\ and \\P\/S] - \\P\J *\ in this and the 
following sections, we denote these by Np and Mp respectively. 

The basic idea of our algorithm is to use the prune and search strategy of Megiddo [21]. 
Clearly, we cannot hope to compute CEITRE(P) by a naive application of this 
technique, since the centre of a reduced set need not be the same as the centre 
of the original set. However, it might be possible to prune points in such a way 
that the centre of the pruned set is a subset of the centre of the original set. If so, 
by repeated pruning we may at least be able to compute a centrepoint, if not some 
larger subset. Below we show that this is indeed possible, and as a first step towards 
this goal we make the following important observation. 
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Figure 23: CEITRE(P - T) C CEITRE(P) 

Observation 4.1 IfTis the set of vertices of a triangle that contains 
then CEITRECP) is a subset o/CEITRE(T U P). 



Proof: Let c be any centrepoint of P, i.e. c 6 CEITRE("P). By definition, any 
closed halfplane, say //, that contains c also contains at least N-p points of P. Since 
by assumption c is contained in the triangle formed by T, we have T O H -^ so 
that H contains at least one point of T '. Thus H contains at least Nj> + 1 = 
points of PUT. Since // is arbitrary, it follows that c is in CEITRE(P U T). 



The above observation has the following important consequence. If we can find a 
set of three points, T, in P such that the triangle formed by these points contains 
the centre of P T, then by discarding these three points we can achieve the goal of 
ensuring that the centre of the pruned set is a subset of the centre of the original set. 
The following lemma gives a sufficient characterisation of such a triplet of points. 

Lemma 4.1 Let T be three points ofP such that HULL(7Vp - l,P) is contained in 
the (closed) triangle formed by T. Then CEITRE(P - T) is a subset o/CE!TRE('P). 



Proof: Let c be a centrepoint of P - T and T be the triangle formed by T. 
We claim that c lies inside T. Otherwise, if c lies outside 7\ and therefore outside 
VULL(Np 1,"P), then there exists an open halfplane that contains c and at the same 
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r' 
Figure 24: Removal of three points may expand the centre 

time contains less than Np - 1 points of P. It can be easily seen that this halfplane 
contains less than Np - 1 = Np. T points of P - T (Fig. 23). This contradicts the 
assumption that c is a centrepoint of P- T. Hence c is contained in T and therefore 
CEITRE(P T) is also contained in T. 

The proof of the result now follows from Observation i.l. | 



Remark. We would like to point out a subtlety involved here. Had we chosen 
the triangle T to contain CEITRE('P) instead, we could not have guaranteed the 
conclusion of the above lemma. Fig. 24 shows why. 

The above lemma suggests an algorithmic solution to the problem of computing a 
triplet of points that can be pruned. Since an open halfplane that contains less 
than k points of P does not interesct HULL(fc, P}, we find three open halfplanes, each 
containing less than Np 1 points of P and situated so that the intersection of their 
complements is a bounded triangle. This triangle contains EULL(Np - l,P). If this 
triangle is of non-zero area then a required triplet is formed by choosing a point 
each from the closure of the pairwise interesctions of these halfplanes (Fig. 25). 

The snag in this solution is that there are configurations of points for which we 
cannot find such a triplet for any choice of these open halfplanes. An example of 
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Figure 25: Pruning of triplets T from P 
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Figure 26: A pathological configuration 

such a configuration is shown in Fig. 26, where the points are evenly arranged on 
the circumference of a circle. 

To overcome this problem, we enlarge the scope of the above lemma, allowing for the 
choice of four points. For this we briefly review the concept of a Radon point. Any 
set of at least four points in the plane can be partitioned into two disjoint subsets 
such that the intersection of their convex hulls is non-empty. A Radon point of this 
set is a point in this intersection. A Radon point of four points is unique when these 
points are vertices of a quadrilateral of non-zero area (Fig 27). 

Lemma 4.2 Let Q be any four points ofP such that the (closed) convex hull of Q 
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Figure 27: Radon point(s) of four points p,g,r and s 

contains HULL(JV> - l,P). Then CEITRE((P - Q) U {q}) is a subset o/CEITRE(P), 
where q is the Radon point of Q. 

Proof: Let c be a centrepoint of (P - Q) U {q}. Consider any closed halfplane 
H that contains c. Then, by definition, it contains at least N(p-Q)u{i} points of 



We claim that c lies in the convex hull, Q, of 2. Let p be a point that lies outside Q, 
and therefore outside E\JLL(Np 1,P). It is then possible to find an open halfplane 
that contains p and contains less than Np 1 = A^(p.g)u{,} points of (P Q) U {q}. 
Hence p is not a centrepoint of (P - Q) U {q}. Therefore c cannot lie outside Q. 

To complete the proof, we have to show that any closed halfplane H which contains 
c contains at least Np points of P. Clearly, H contains at least Np-Q)^ points 
of (P - Q) U {q}. Three different cases arise, depending on the relative positions of 
the points of Q. 



Case 1: The four points in Q form a non-convex quadrilateral. 

This case is a trivial application of Lemma 4.1. The three convex vertices 
of Q form a triangle that encloses mLL(Np - 1,P) and the concave 
vertex is q. Thus by Lemma 4.1 CEITRE((P - C) U {q}) C CEITRE(P). 
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Figure 28: Substitution of Q by its Radon point q 

Case 2: The four points of Q form a convex quadrilateral but their Radon point 
q does not belong to H . 

Since the quadrilateral Q and the halfplane H both contain c, their 
intersection is non-empty. Thus H contains at least one of the vertices 
of Q and therefore at least Np = -/V(p_Q)u{ f } + 1 points of P. 

Case 3: The four points of Q form a convex quadrilateral and their Radon point 
q belongs to H (Fig. 28). 

In this case H contains at least two points of Q. We can therefore delete 
q from H and still claim that H contains at least N-p points of P. | 



Thus in all cases H contains at least Np points of P. Since H is arbitrary, c is 
a centrepoint of P as well. Hence CEITRE((P - Q) U {q}) is a subset of CEITRE(P). 



The above lemma is the cornerstone of our pruning mechanism. In the next section 
we will show how to use ham-sandwich cuts to make a clever choice of four open 
halfplanes so that we can prune a fraction of the input set by repeatedly applying 
the last two lemmas. 
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4.3 How to prune 

In this section and afterwards, we use the words left, right, up and down, wherever 
these are unambiguous, to simplify the arguments. 

Suppose we choose four open halfplanes, call them L, U, R and >, (mnemonics for 
Left, Up, Right and Down respectively) such that each contains less that Np - 1 
points of V and its closure at least Np points, and together they enclose a non-zero 
bounded area. Why do we expect this choice to give us a triplet/quadruple of points 
Q satisfying the conditions of Lemma 4.1/4.2 ? We give an intuitive justification of 
this below. 

If the pairwise intersctions of "adjacent" halfplanes ( i.e. L and (/, U and R, etc.) 
are empty we would get a configuration as shown in Fig. 29. In this configuration, 
the total number of points in all the halfplanes taken together exceeds the total 
number of points in P by approximately one-third! This is impossible. So we 
might attempt to construct the four halfplanes in such a way that this excess 
is distributed evenly among the pairwise intersections of the adjacent halfplanes 
and thereby obtain approximately Mp triplets/quadruples of points satisfying the 
conditions of Lemma 4.1/4.2. 

It is possible to do this as the construction below shows. 



4.3.1 Computation of Open Halfplane L 

We fix L as follows. We determine an extreme point p of P with minimum abscissa, 
and join all the remaining points to it. We compute the line that passes through 
p such that its slope is the Np - 1-th largest of the slopes of the above lines. The 
open halfplane above this line is chosen to be L. Clearly, it takes linear time to 
compute L. This way we make sure that L contains less than Np - 1 points and its 
closure contains at least these many points of P. Moreover, this halfplane contains 
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Figure 29: The intuition behind pruning 

at least two points on its boundary. Really speaking, this latter requirement is not 
necessary but it helps us to treat the four halfplanes uniformly in the analysis of the 
algorithm. 



4.3.2 Computation of Open Halfplane U 



Since the point set P can be degenerate we need to be careful in the construction 
of U in order that none of the closed quadrants determined by the boundaries of L 
and U contains too few points of P. This will also result in an even distribution of 
points among the pairwise intersections of adjacent halfplanes. 

To achieve this we use the ham-sandwich cut algorithm of Megiddo [23]. However, 
the ratios in which we propose to divide the point sets are arbitrary. As we show 
below, Megiddo's algorithm can be easily adapted to take care of this aspect. 

As is usual, we consider the dual problem, letting the bondary of L be the y-axis in 
the primal plane. Under the duality transformation that we consider, points that 
lie on the boundary of L map to horizontal lines; the set of points that lie in L 
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map to a set of lines, _, with negative slopes; and the set of points that lie in 
interior of L map to a set of lines, + with positive slopes. In order to get the type 
of ham-sandwich cut we want, we put the horizontal lines in _. Let us assume 
that + contains at least Np Mp lines in the dual plane. It can be seen that the 
boundary of U corresponds to that point in the dual plane which has less than Mp 
lines of C~ above it; at least these many lines of C~ passing through or above it; 
less than Np - Mp lines of + above it; and at least these many lines of + passing 
through or above it. We compute this point as follows. 

It can be easily seen that we can resolve a query for any positive (negative) slope 
query line in Megiddo's method as follows. We first compute the Mp-ih (Np - Mp- 
th) intersection of the lines in C~ (+) with the query line. We then count the 
number of lines in + (C~) lying strictly above and the number of lines in + (~) 
passing through this intersection point. If the sum of these two counts is smaller 
than Np - Mp (Mp) then the solution point lies below the query line. If the first 
count is greater than or equal to Np-Mp (Mp) then it is above the query line, else it 
is same as the intersection point. By thus changing the method of query resolution, 
we get the solution point in linear time. This gives us the required halfplane U in 
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the primal plane (Fig. 30). 



4.3.3 Computation of Open Halfplanes D and R 



The open halfplane D is determined with respect to L in a similar manner by 
changing "above" to "below" throughout the above discussion. Thus we make sure 
that D contains the "lower end" of the boundary of L whereas U contains the "upper 
end". 

The halfplane R is also determined similarly except that U plays the role of L here. 
We also ensure that R contains the "right end" of the boundary of U whereas L 
contains the "left end". 

The idea behind choosing D and R in this way is to make the boundaries of the 
haifplanes U ', R, D and L the adjacent sides of a bounded quadrilateral such that 
the haifplanes face "outwards". This is so because if the interior of ~L contains at 
least Np Mp points of P then the boundaries of both U and D are disjoint from 
the boundary of L and it can be seen that the intersection of the complements of 
the above haifplanes is bounded. As explained before, it can also be seen that RC\ D 
contains non-zero points of P. Moreover, since each of these half planes contain 
less than Np 1 points of P, the intersection of the complements of these encloses 
HULL(7V>-1,:P). 

However, in the computation of U it is quite possible that the set + contains less 
than Np - Mp lines. The consequence of this is that the computed U may have 
the same boundary as that of L and thus the intersection of the complements of the 
open half planes may be unbounded. This is inadmissible in our algorithm. So we 
need to take care of this degenerate case separately. 
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4.3.4 The Degenerate Case 

If there are more than Np points on the boundary of L then the boundaries of U 
and >, as computed above, is the same as that of L. This degeneracy is detected in 
the algorithm when + contains less than Np - Mp lines i.e. interior of L contains 
less that Np Mp points of P. There is no loss of generality if we assume that L 
also contains less than Np Mp. Otherwise, we can switch the sides of L and let the 
interior of L be our new L. We will then be able to compute the other halfplanes, 
as required, with respect to this open half plane. 

Let the interior of T be the open half plane R. We compute the open halfplanes V 
and D as follows. We first distribute the alternate points of P on the boundary of 
L between the sets R and L. Let these resulting sets be SR and S L respectively. We 
then compute U and D by the ham-sandwich cut algorithm such that these satisfy 
the following properties: these contain less than total Np - 1 points of P and contain 
an equal number of points of S R and S L . The halfplanes U and D are computed such 
that these contain the "opposite ends" of the boundary of L. It can be seen that if 
L, R, U and D are computed in this manner then each of the pairwise intersections 
of the adjacent halfplanes contains at least Mp (approximately 2Mp) points of P. 
Since the intersection of the complements of these contains WU^Np - Mp,P], we 
can apply pruning at least Mp times similarly as in the non-degenerate case that 
we describe below. 

The degeneracy of open halfplanes taken care of, we may assume safely that the open 
halfplanes I, U, R and D can be computed such that these meet our requirements. 



4.3.5 The Pruning Step 

We have been able to ensure by the construction of the halfplanes as above that the 
closure of each of the sets L 17, L D and R U contains at least MP pomts of 
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P. We shall prove later that the closure of R fl D also contains at least these many 
points. However, for the rest of this section we will assume this. 

It is now clear how we can prune points. Two points of detail, however, must be 
noted. First, in order to ensure that the conditions of Lemmas 4.1 and 4.2 remain 
valid throughout the pruning step, we must choose a triple or a quadruple of points 
in such a way that, whenever there is a possibility that the conditions of the above 
Lemmas are violated in the successive pruning steps, we delete an interior point in 
an open halfplane. Second, to maximise the number of points that are pruned, we 
must ensure that no two points of a quadruple, selected for pruning, belong to either 



To implement the above observations we maintain the points that are candidates for 
pruning in six disjoint sets, viz., PLU, PVR, PRD, PLD, PLR and PUD- The points 
on the boundaries are put in the relevant sets. So, the four sets, LC\U, LC\D, RC\U 
and R C\ Z), are now effectively divided into six sets, three of which correspond to 
choices of triangles T. 

We discard the triangles T and substitute quadruples Q by their Radon points 
such that a maximum number of the above half planes contain an interior point. 
Substitution of Q is done as follows. If Q form a convex quadrilateral we delete it 
from P and add the intersection point of the diagonals to P. Otherwise we delete the 
convex vertices but retain the concave one. We can repeat this pruning procedure on 
the reduced set of points thus obtained, since the halfplanes L, U, D and R continue 
to contain less than Np - 1 points of the reduced set P, until one of the four sets is 
empty. We note that this reduces the size of P by approximately one fourth. 



4.4 The Centrepoint Algorithm 

It is now clear from the discussion in the previous sections, how we can find a 
centrepoint of P. 
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In each iteration we compute the points that are to be discarded or replaced. By 
throwing away these points we reduce the size of the set by a non-zero fraction. 
When the size of the set becomes so small (of size at most 10) that no more points 
can be discarded we halt the pruning procedure and compute a centrepoint by any 
straight forward method. 

The algorithm is given below. 

Algorithm 4 CENTREPOINT 

Input: Set of Points P 
begin 
do 

Compute the open half planes L, U, D and R 
Update P by deleting T and replacing Q by 

their respective Radon points 
while there is some replacement /deletion 
enddo 

Compute a centrepoint by any bruteforce method 
end | 

We justify that anytime during the pruning step each halfplane contains less than 
NK 1 points where K is the current set of points. 

We first argue for a deleted point that is also an interior point of a halfplane. We 
first consider the case when four points are pruned and their Radon point is added 
to the set P. If a point that is pruned lies in the interior of a halfplane and the 
Radon point does not lie in that halfplane then the number of points in this half 
plane is decreased by one. Since the total number of points decreases by three, the 
number of points is less than N*-l=N T -2 where ft is the new set. Now suppose 
that the Radon point also lies in this halfplane. Since a halfplane that contains a 
Radon point of four points contains at least two of these points, in this case also the 
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Figure 31. The worst case of pruning! 

number of points in the halfplane decreases by one. The same argument holds for 
triplets of points. 

Now we consider the case when there is no interior point among the four pruned 
points in a halfplane. In this case the Radon point also does not belong to it. We 
consider the worst case in which for every choice of a quadruple or a triplet of points 
those chosen from this halfplane lie on its boundary. This case is simpler to analyse 
and it is easy to see that same argument extends to the other cases. It is clear that 
when we start pruning, the number of points in each of the open halfplanes is less 
than Np 1 points. Let us see what happens when we have to prune the last triplet 
or quadruple of points after applying pruning Mp 1 times. The maximum number 
of points in this open halfplane is less than 

Np-Mp 

since Mp 1 points in the corner quadrants are already pruned. 
The total number of points at this instant is 

\P\ - ZMp + 3 



The number in the previous expression is exactly one less than the ceiling of this 
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number. Thus the last set of points can also be pruned. It can easily be seen that 
in the intermediate steps also we can do the pruning. 

The proof of correctness of this algorithm now follows from Lemmas 4.1 and 4.2. In 
the next section we give an analysis of the running time of this algorithm. 



4.5 Analysis of the Centrepoint Algorithm 

For the purpose of the proofs in this section, we assume that the boundaries of the 
halfplanes L,U,D and R do not contain any point of P. This assumption is not 
necessary to establish the linearity of our algorithm but it simplifies the proof to a 
great extent. This can be achieved by slightly perturbing the points of P such that: 
no point migrates across the boundary of any halfplane; each corner region, such as 
LOU etc., contains exactly Mp points of P\ and each halfplane contains at least 
Np points of P. Intuitively, such a perturbation does not matter because if we 
can prove that the perturbed set PRO contains Mp points then these many points 
also belong to the closure of PRD in the non- perturbed set P. As a consequence 
we do not need to state explicitly whether the halfplanes are open or closed in the 
following discussion. 

To prove that the algorithm is linear we have to show that the size of P is reduced by 
at least a fraction in each pruning step. We know by construction that each of the 
sets, PUJ, PUR and P LD , contains Mp points. We have to show that PRD also con- 
tains at least Mp points in order to prune at least these many triangles/quadruples 
from P. For this we will have to consider all the possible relative positions of the 
halfplanes L, U, R and D. Several cases arise. A few of these can be straight away 
discarded by using the following lemmas. 

Lemma 4.3 Let F, G and H be three halfplanes. Then F H 5 n # is a bounded 
triangle if and only if F fl G C H. 
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Figure 32: Intersection of three halfplanes 
Proof: Straight forward(Fig. 32). | 

Corollary 4.1 T fl <5 fl H is a bounded triangle if and only if F G C H if and 



The sets t/, , R and D satisfy some additional constraints also on account of their 
specific methods of construction. 

Let us denote by PQH the intersection point of the boundaries of any two halfplanes 
G and H. For the sake of simplicity, we may assume without any loss of generality 
that no three boundaries of the above halfplanes intersect at a point. If they do, 
then we can treat the said configuration in one of the cases discussed later on. 

Lemma 4.4 The intersection of the halfplanes U and D is contained in L if and 
only if PUD w contained in L, i.e. U fl D C L <$=> p UD G L. 



Similarly, L fl R C U <=* pm U . 

Proof: We prove only the first equivalence. The proof of the second is similar. 
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(=>) Easy. 

(=) Let PUD be in L. Then four cases arise, depending on the orientation of the 
half planes U and D. 



Case 1: 

Case 2: U r\ D C L 

These cases are not possible. For by construction, the halfplanes U and 
D contain the "opposite ends" of the boundary of L. 

Case 3: U O^ C L 

This case is also impossible since it implies that the centre is empty. 

Case 4: U D D C L 

This is the only permissible case(Fig. 33). I 

Hence proved. I 

We have the following similar lemma when pvo lies in L. 

Lemma 4.5 The intersection of the halfplanes U and D u contained in I if and 
only if PUD w contained in T, i.e. U D C L <= PUD L. 



Similarly, LHRcU <=> PLR U. 

Proof: Here too we prove only the first equivalence. 

(=>) Easy. 

(=) Let PUD be in I. Then four different cases arise. 
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Figure 33: Either UnDcLoiUnDcL 

Case 1: U nT) C~L 

Case 2: Ur\D d 

These cases are not allowed by construction. 

Case 3: UnDcL 

To see that this case is also not permissible, let us count the number of 
points in P L . First we prove that P L - P LU U P LD . 



PL = 

= PLUDUP LU D U ' P LUD' since Pjffjy = 0, 



since PWD U P LUD = "PLD and VUJD U P Lt /iy = PLC/- Hence, 

\PL\ = 



U\ + \PLD\-\PLUD\ 
-l- A/^ - \P LVD \ 



79 



This contradicts the fact that P L contains at least N P points. Thus this 
case is also not possible. 

Case 4: U O D C I 

This is the only permissible case. | 

Hence proved(Fig. 33). | 

The above lemmas have the consequence that either U D C L or U C\ D C L, and 
similarly either LClRcUoTLnRcU. Now we can prove the following theorem. 

Theorem 4.6 There are at least Mp points in PRO- 

Proof: For the proof, we again consider all the possible relative positions of the 
halfplanes L, /, R and D. The following cases arise, depending on which of the 
four quadrants formed by the boundaries of U and L contains p/u). 



Case 1: PRO e U D L (Fig. 34) 

In this case D D U H L and R U fl L are non empty. Therefore, from 
Lemma 4.4, 

UC\D CL and 



Using Corollary 4.1 these respectively imply that 

and LC\DcU. 



Let x be a point of R fl D. Then x lies in exactly one of the sets, UftL, 
C/nI,FnlorFnI. Since UC\RC 1 therefore x <?U C\ L. Similarly 
x # U H I. Now R H D is convex and there is a point p RD that lies in 
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Figure 34: p RD (7 H L 

U fl L, therefore x does not lie in U L either because a convex region 
cannot intersect only the opposite quadrants of a pair of lines. Thus 

RC\DC UC\L, 
and hence, by Corollary 4.1, 

Ur\RC~D and 

Now, 

V = 



= PLUPi u V'PuR (J ' p lD [J 'PlU'K'B since In /? C I/ and InZ? C 
As these are disjoint sets 
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V 





Figure 35: PRO U fl L 



> in 

which is a contradiction. Hence this case does not occur. 

Case 2: PRD TJ n I (Fig. 35) 

In this case R n L (7 * 9 therefore RC\LC~U (Lemma 4.5). We will 
have to consider four subcases. 



1. RKDC L 

2. JROTJ C I 

Since RKDCL therefore fl and Contain the "opposite ends" of 
the boundary of L. Since RHLcV therefore from Corollary 4.1 
R and C/ contain the "opposite ends" of the boundary of L. This 
implies that D and U contain the "same ends" of the boundary of 
L, which is not possible by the construction of D. 

3. ~RCiDC L 
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By Lemma 4.3 this implies, 

~L O D C R. 
Therefore, 



Thus 



> N T - M P . 

4. RnT> C L 
This implies 

RC\U C D. 

Otherwise consider a point x 6 R Pi U O D, 



Again, 

xe/2 => x e ROL 

which is a contradiction. 
Therefore, 



Hence, 

\PRD\ > \Pmj\ 
> M P . 

Since PRO L so these four cases exhaust all the possibilities. 
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Figure 36: p RD U n L 

Case 3: PRD 6 U fl I 

This case is similar to the previous one. 

Case 4: p RD 17 n I (Fig. 36) 

Again four different subcases are possible. 



1. RrdcV <uidUr\DcL 

Consider a point xRr\D. lix L then x U by fln I C U. If 
x G I then xeI7by[/ncL. Since R n D is convex, it cannot 
intersect exactly two opposite quadrants of a pair of lines. As there 
is a point PRD in U H I hence R H D C U 1. 
Assume that |P/jp| = m. Then, 



since /? n D c I/ 

As the sets on the left of the last equality are disjoint sets, 
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Substituting the values of the different terms, we get, 



< m + |P| - 4A^ 4- 3M, 



Since the size of the set Pyj^ is non-negative, 
m > Afp. 



2. RDLCU an 

Since PRO 17 H T therefore PRD lies in U. There are four different 
possibilities. 

t RC\DCU 

RKDcU 

By a similar argument as in case 2 we can show that the above 
two cases are not possible. 

RODCU 

Consider a point x R fl L. 



U => xe t/ 

Rr\~DcU = x 



=^ x 6 I n I, 

which is a contradiction. Hence flfll = 0. Therefore #01 C 

and U C\ D C L a subcase we shall consider later. 

RC\DCU 

Let |T>UD| = m. We compute | 

Now, 

P = 
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= P D U PO R U PW U PDRUL U ^nra, since TEH 7 C 27 



since L D / C #. 

Substituting the values of different terms in the above equation 

we get 

\P\ > Np + Np-m + Np-Mp + Np 



Since the size of the set P"D^UX ls non negative therefore, 
m > Mp. 

3. ^OLcZ/and^nDcL 

This case can be dealt with in a similar way as above. 

4. RC\LCU andUr\D cl 

Since PUD / here also we shall consider four different cases. 

RC\D CU 

RC\DCU 

These cases do not occur for the same reasons as discussed 
above. 

RHDCU 

From Corollary 4.1, 

RC\U C Rf\D. 
Since R O U contains Mp points of P, 

> \P*u\ 

> Mp 
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RnDcU 

We need only consider the case R D D C T since other cases, 
viz. Rr\D C~l, SflDcIandBnUcI, can be dealt 
symmetrically as above (actually they leads to a contradiction 
when considered along with R C\ D C U), Let PRD contain m 
points. As before we compute the number of points in 



= ?D U Pp,, U Pjfr U PBKPJ, U PZJRPT, since fin UC D 



since I D D C fl. 

Substituting the values of different terms in the above equation 

we get 

\P\ > Nr + NT-m + Nr-Mr + Np 



Since the size of the set Pjji'R'D * s non negative therefore, 

m > Mp. 
Thus theorem is proved for the last case PRO Z7 fl L also. | 

Hence proved. | 

Combining the earlier theorems and lemmas we get the following result. 
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Theorem 4.7 A point in the centre of a set can be computed in linear time. 



Proof: In each iteration at least 3Mp (~ |P|/4) points are deleted. 

If T(n) is the running time of the algorithm for an input set of size n (\P\ = n), 
then it satisfies the following recurrence. 

T(n) < 
=* T(n) < 



Since T(n) = 0(n) from the above recurrence relation, the claim of the theorem 
follows. | 



4.6 Concluding Remarks 

We have presented an optimal algorithm for computing a centrepoint of a finite set 
of points in the plane, thus providing one more example of the power and versatility 
of the prune and search paradigm. 

It would be worth exploring how this speeds up algorithms which uses the centre- 
point computation as a basic subroutine. 



Chapter 5 

Designing Algorithms Using 
Partial Sorting Networks 



In this chapter we present a general technique that is a sequel to the parametric 
searching technique of Megiddo. The latter technique is used to design serial 
algorithms with the help of efficient parallel ones. However, these algorithms are 
not very efficient for the class of problems that use parallel sorting networks in their 
solutions. We modify this technique so that we obtain optimal linear time algorithms 
for some of these problems. We do this by synthesising prune and search technique 
and parametric search technique. 



5.1 Introduction 

The parametric searching of Megiddo is as follows Let there be an efficient parallel 
algorithm for a problem A such that solution of A can be used in the solution of 
another problem B. Then, in some cases, we get an efficient sequential algorithm 
for B by exploiting the efficient parallel mechanism of the parallel algorithm for 
A. This technique has been applied to a wide variety of problems yielding efficient 
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algorithms. In particular, we achieve good results for the parameterised problems 
that use parallel sorting algorithms in their solutions. We do this by replacing 
the evaluation of the parallel comparisons of an iteration in the parallel version by 
simultaneous resolution of these in the serial version. The running time of these 
algorithms is further improved by introduction of weights in the comparisons [8]. 
We then simultaneously evaluate at least a fraction of total weight of these in every 
iteration to design more efficient algorithms. 

We further improve the running time of the above algorithms by introducing, wher- 
ever applicable, prune and search in these. We do this by seeking to compute the Jb-th 
largest element of the input set instead of seeking to sort it as in previous techniques. 
For this, we run the sorting algorithm on a given input for a few iterations and then 
prune the set. It can be easily seen that this approach is useful only where pruning 
is applicable. 

This chapter is organised as follows. In section 5.2 we review a few definitions and 
concepts related to the AKS sorting networks. In section 5.3 we present a method 
of applying prune and search using these networks. In section 5.4 we optimally solve 
the problem that was first posed by Megiddo. In subsequent sections we apply this 
technique to linear programming problem, computation of ham sandwich cuts and 
computation of centre points in d- dimensions. 



5.2 Preliminaries and Definitions 

The AKS sorting network sorts a sequence of reals in O(nlogn) comparisons such 
that in each parallel iteration 0(n) comparisons can be performed so that there 
are in total O(logn) iterations. A comparison operation in the network consists of 
comparing the contents of two registers and interchanging these if their specified 
order does not match. 

We can view the process of sorting in the AKS sorting network as the movement of 
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registers in a complete binary tree of n leaves. The order of the leaf (topmost) nodes 
of this tree corresponds to the ordered sequence of n numbers. All the registers are 
initially in the root (bottommost) node of the tree and in every iteration, most of 
the registers in a node move up whereas only a small number of registers move down. 
When each of the registers reaches a leaf such that each leaf has only one register 
present in it then the computation is said to be over and the elements in the leaves 
are in sorted order. 

We use the same notation as in [1] by Ajtai et al. Let K be a set of n registers, 
which contain elements of the input set, in the AKS sorting network. Let the set 
of registers assigned to a node t of the complete binary tree in a-th iteration be 
denoted by S Q (t). Thus, S(t) is H if t is the root of tree, otherwise it is (empty 
set). The restriction of the function S a to a fixed level of height h is denoted by 
S' H . S is defined in such a way that each of the nodes at any height, /i, contain 
same number of registers. Let N(S a ' h ) be this number. We bound this number in 
the following theorem. We reassign the constants 71, q? and c\ used in the paper [1] 
as g, Q and c respectively in the following discussion (Q -C q -C 1/c <C 1). 

Theorem 5.1 ([!]) The number of registers in any node at height h in a-th itera- 
tion has the following properties: 

1. N(S Q ' h ) < q- h 2- a n. 

2. if q- kl 2- Q - l n < c then N(S' h ) = otherwise N(S Q ' k ) > q- h 2- Q - l n. 



We can see from the above theorem that the number of registers in any node in 
a-th iteration is bounded, both above and below, by terms that are in decreasing 
geometric sequence. The AKS sorting network has another useful property. The 
ordering of the elements after each iteration is an approximation of the sorted order 
of the input set. Moreover, at any instant the number of the elements having a 
large displacement is very small. We bound this number in the following theorem. 
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Let the relative position of a number x in a sequence be defined as the fraction of 
numbers in the sequence that are smaller than z. We also assume that the positions 
to right correspond to the larger elements. 

Theorem 5.2 ([!]) Let G^t), where t is a node at height h (h < a), be the set 
of registers in the nodes, right of t, at height h such that the relative positions of 
elements in these are smaller by at least p than the relative position oft. Similarly, 
let H^t) be the corresponding set of registers with larger elements on the left oft 
such that the relative positions of the elements contained in them are larger by fji 
than the relative position oft. Then, 



and, similarly, 

where M = N(S a ' h ) andm>\. 

Now, we seek to apply the Megiddo's technique to optimisation problems with the aid 
of these theorems. In these, we frequently need to prune redundant hyperplanes in 
E d . The multi dimensional prune and search procedure proposed by Megiddo is used 
for this purpose [22]. According to this, we can prune a fraction of the hyperplanes 
by locating the optimal point with respect to a finite set of query hyperplanes. The 
following theorem establishes the existence of these query hyperplanes. 

Theorem 5.3 ([22]) For any dimension d there exist constants A = A(d) and 
B B(d), with < B < l /i, such that A queries suffices to determine the position 
of optimal point relative to at least Bn of n hyperplanes in E d , where 



and 

B(d) = 2 1 - 
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It is clear that this procedure of pruning hyperplanes is applicable only if there 
exists an oracle, say fi, that takes linear time to determine the relative position of 
the optimal point with respect to a query hyperplane. 

In the next section we discuss the pruning of elements in the input set. 



5.3 Pruning Using the AKS Sorting Networks 

Let $ is a set of n elements. Then, there are altogether TV (= ^2) distinct 
comparisons possible between pairs of the elements of S. We assume that the 
output of any comparison depends on the value of a vector x i.e. it depends on 
the membership of x among disjoint subsets of the domain. In the problems we 
discuss, these are the regions determined by hyperplanes. We call these critical 
hyperplanes. If the optimal points x* is known then the complete order of elements 
of S is known. The critical hyperplane of each unresolved comparison intersects the 
localised region. Furthermore, it is to be noted that localisation of x* with respect 
to one of these hyperplanes as well as localisation with respect to O(n) of these 
(only a fraction of these will be resolved at a time using multi dimensional prune 
and search) require linear time( Theorem 5.3). 

We pose a hypothetical problem of computing the Ar-th largest element, s*, of 5 at 
x*. What are the elements that can be pruned ? Clearly, these are the elements 
that are not the fc-th largest element at x*. We formulate our problem as follows: 

Subproblem I How can we determine, in linear time, O(n) elements ofS that are 
not same as s^ given an oracle Cl that takes O(n) time to answer the following query 
which side of a given hyperplane does x* lies? 

The real crux of the problem now is to identify the elements that have large dis- 
placement from a*- For determination of these, we generalise the method of Lo and 
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Steiger in [18]. We determine a prismoid, whose base is the localised region, about 
the *-th largest element such that most of the elements do not intersect it. The 
elements that do not intersect this prismoid are always at a non zero finite relative 
displacement from s k in the localised region and therefore, can be pruned. 



5.3.1 What to Prune ? 

Let S be a set of hyperplanes in E d . Suppose we can localise x* in a region J such 
that the total number of unresolved comparisons is only a small fraction of N. We 
construct an infinite prism with J as its base such that only few of the intersections 
of any pair of elements of 5" intersect J. Then we prune the hyperplanes by the 
following theorem. 

Theorem 5.4 Let V be the prismoid determined by k - (d - \)\en\ and k + (d - 
l)fen] largest intercept on each vertex of localised region (which is a simplex). If 
there are less than "t^fenl 2 unresolved comparisons in the localised region then we 
can prune approximately (n-d[^en\) elements in the computation of k-th largest 
element. 



Proof: Each unresolved comparison corresponds to the intersection of a pair of 
hyperplanes in 5 which also intersects the infinite prism of which P is a part. 

We compute the minimum number of these intersections needed so that the surface 
formed by fc-th largest elements (fc-level) may lie above or below the prismoid at an 
interior point p. We only analyse the case when the surface lies above P at p. 

Let us consider the following worst case configuration. In this configuration, each of 
the hyperplanes that is above the fc-level and below the k + (d - 1) \en\ -level at any 
of the vertices u t , with < t < d of the localised region is above the prismoid at p. 
Let the set of these hyperplanes be H. Also, in the minimal case, no hyperplane that 
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is below the Jb-level at any of the t\'s is above the prismoid at any point. Otherwise, 
the number of intersections in J increases. Similarly, no hyperplane, that is not in 
H but which is above the k + (d 1) [en] -level at any of the t;,'s, intersects prismoid. 
Furthermore, none of these hyperplanes intersects any other hyperplane over J. 

Now, the hyperplanes in H in this configuration can be divided into several equiv- 
alence classes defined as follows. If two hyperplanes intersect or do not intersect 
the prismoid at the same vertices then these belong to same equivalence class. 
Each equivalence can be represented by a set of vertices of the base, at which the 
hyperplanes of this class intersect P. We call these as intersecting vertices of an 
equivalence class. A hyperplane in one equivalence class intersects a hyperplane of 
another equivalence class if their intersecting vertices are not related by the relations 
of either subset or superset. 

If there are d equivalence classes of size \en] such that each has only one intersecting 
vertex then the number of intersections above J is minimum. The count the number 
of intersections in this configuration is "t^F^] 2 . If tne number of intersections is 
less than this then the fc-level will lie inside the prismoid. 

We can prune those hyperplanes which are above or below the prismoid at all the 
vertices v,. The maximum number of hyperplanes that will intersect the prismoid is 
when the equivalence classes are chosen as above and the intersections are divided 
equally between the hyperplanes both above and below the Jt-level. Hence, we can 
prune the rest of approximately (n d [v/2enj) hyperplanes in 5. | 

We can prove another similar theorem. 

Theorem 5.5 If there are less than ^fen] 2 unresolved comparisons in the feasible 
region then the prismoid is formed by the largest (smallest) and (d l)[en] largest 
(smallest) intercept on each vertex of feasible region (which is a simplex) and we can 
prune (n d\en\) elements in the computation of the largest (smallest) element. 
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To use prismoid method we need to ensure that number of unresolved comparisons in 
the feasible region should only be a small fraction of N. Then we can guaruntee the 
existence of a non zero fraction of prunable elements. In next section, we compute 
the minimum number of iterations of AKS sorting network needed to ensure this. 



5.3.2 Why can we Prune ? 

We use the properties of AKS sorting network stated in the previous section to 
prove that the number of unresolved comparisons in a-th iteration is bounded by a 
number that decreases in a geometrical progression with a. 

Theorem 5.6 The number of unresolved comparisons in the AKS sorting network 
in the a-th iteration, out of a maximum of"Ci ones, is bounded from above by 



Proof: Let us first count the number of unresolved comparisons of an element x 
in a register present in a node at height h in a-th iteration. Let p be 2/c and t be 
that node which has the same relative position as x. Then 
2 m /M < ji, 
2 m < pM 
< 2M/c 



Since fl*-*2 ! n > c, therefore, m = 1 satisfies the above inequality. 

Thus, the number of unresolved comparisons of elements at height h with a larger 
relative displacement than p from x is bounded by 



This quantity when summed over each level of the computation tree is at most 



The number of remaining unresolved comparisons of i with elements having relative 
displacement smaller than \i is at most 2/i = 4n/c. We take half of this in the total 
as each of these unresolved comparisons is counted twice in the sum. Further, we 
also add unresolved comparisons corresponding to any two registers in the same 
node to the computed total. Therefore, number of total unresolved comparisons is 
at most 

)2~ a nQ -I- 2n/c ) n + ] 5 " 

~~ ? / Q<h<a 



Hence proved. 

We can make this fraction as small as we need by choosing a, g, (J, c appropriately. 
We can thus bound the number of unresolved comparisons in each iteration. But 
there still remains an algorithmic problem. Resolving all comparisons of even 
only a first few iterations of the AKS sorting network requires O(logn) steps of 
Megiddo's multi dimensional prune and search. This will take at least O(nlogn) 
time. However, this is not agreeable. So, in order to make the pruning procedure 
linear, we have to leave another small fraction, say 6/2, of comparisons unresolved. 
We observe that this adds at most 2(6n/2)n in the above quantity i.e. two elements 
for each unresolved comparison (and these may remain uncompared with all the 
other elements). Then, the total number of unresolved comparisons is 



2-q* 
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where 6 is also chosen appropriately. Now, we can apply Megiddo's multi dimen- 
sional prune and search on the critical hyperplanes until only 6n/2 comparisons of 
the total comparisons of first a iterations are left unresolved. Clearly, this takes 
linear time. 

The method of application of the technique will become more clear in section 5.4 
when we apply our technique to a few problems. 



5.3.3 How to Prune ? 

We can view AKS sorting network as having O(logn) iterations with at most n/2 
comparators in every stage. We give a weight of 4~ J to each comparator at a depth 
j as in [8]. We choose the values of 7, Q, c, 6, a and appropriately and assign 
weights to comparators of first a stages only. The construction of first a stages of 
AKS sorting network takes linear time. We give details of the algorithm to compute 
x* below. 

Algorithm 5 Prune And Search on AKS Sorting Network 

Input: 5 
begin 
do 

Let C be the set of unresolved comparisons 

and W its total weight 
do 

Solve (d-l)-A(d-l) queries using ft ( Theorem 5.3) 

Resolve at least B(d - 1) comparisons 

Update C and W appropriately 
while |C| > 6ri/2 
enddo 
Prune the hyperplanes in 5 ( Theorem 5.4) 



Update k and S accordingly 
while 5 is modified 
enddo 

Apply any brute force method to compute x* 
end | 

In the above algorithm (d - \)A(d - 1) queries are put to the oracle because we 
need to ensure that the base of the prismoid, where x* is located, be a simplex. So, 
whenever a query hyperplane intersect the base, we make at most d 2 additional 
queries to again make it a simplex. 

It can be easily seen that this algorithm runs in linear time if the time complexity 
of oracle is O(\S\). 



5.4 Intersection of Median of Straight Lines with 
a Given Line 

Let there be n non- vertical lines /,, where i = 1, . . . ,n in a Euclidean plane. Each 
of these lines is represented by a linear equation y = /,(x). The median level is 
defined as a function, F, such that F(a) is the ordinate of [n/2J-th intersection of 
the lines with the vertical line x = a. We can evaluate F(x) for any x in linear time 
by computing median of /,(x)'s. 

Let us compute intersection of median level and a given line /, which is represented 
by y = /(x). We are guarunteed a solution if the slope of this line is not the median 
slope if n is odd; does not lie between the [n/2J-th and [n/2 + IJ-th slope if n is 
even. 

For the sake of simplicity we assume that all the lines have positive slope and the 
line/ is the x-axis. Then we have to determine x* such that F(x*) = /(x*) = y* = 0. 
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At x -> -oo and at x -> +00 function F(x) is same as the line of median slope. 
Therefore at x - +00 it is above the x-axis and at x -* -oo it is below x-axis. 
Since function F(x) is continuous, the equation F(x) = has at least one solution. 
This problem was first posed by N. Megiddo to motivate the technique of parametric 
searching. 

In the next subsections we briefly review how Megiddo and Cole applied their 
techniques to solve this problem. 



5.4.1 An Algorithm by Megiddo that Motivates the Tech- 
nique of Parametric Searching 

We first make an important observation. The order of lines changes only at inter- 
section points of a pair of lines. Hence, intersection points play an important role 
in the evaluation of x*. 

Let x l3 be the abscissae of the intersection point of lines /, and /,. 

We can compute x* naively as follows. First we identify intersection points of each 
pair of lines. Then we search for two values x 1 and x 2 such that F(x l ) < < 
F(x 2 ) and there is no other x tj in the open interval (x 1 ,! 2 ). We search for such 
an interval by employing binary search. This requires O(logn) F-evaluations and 
median computation of subsets of x tj whose cardinalities are TV, W/2, #/4, . . . points 
where N = "C 2 - Thus, this algorithm runs in 0(n 2 ) time which is dominated by the 
evaluation of all the intersection points. 

Now we employ Megiddo's technique to design an algorithm that runs in O(n log 2 n). 
We evaluate F(x*) with x* not specified by computing the median of /,'s. The 
median-finding algorithm compares values of function /,(x*) and /,(x*) for different 
i's and j's. The outcome of such a comparison changes only at intersection point 
x . Hence, if F(x ) > then x* < x tj ; if F(x tJ ) < then x* > x,,; and if F(x tJ = 
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then x* = x t -j. Once we know the outcome of a comparison we proceed with the 
evaluation of median. Now the value of x* is restricted in the interval ( oo,x tj ) or 
(xij,+oo). At the end of median-finding algorithm we know the median line /, at 
x = x*. The intersection of this line with x-axis yields x*. Clearly, this algorithm 
is also O(n 2 ). Since the above algorithm is no better then the previous one we try 
another approach using AKS sorting networks. The AKS Sorting network employs 
P = O(n) processors in every stage and sorts its input in O(log n) time. Instead of 
evaluating function F(x) at each of the x,/s, we first wait for all the P values in 
a parallel iteration to be computed. Then we sort these P values and do a binary 
search to find the interval which contains x*. Then, we resolve all the comparisons of 
this parallel iteration. "When all the iterations of AKS sorting network are completed 
then the lines are sorted at x* and computation of x* is easy. 

Since each stage takes O(nlogn) time and there are O(logn) iterations, so the 
overall running time of this algorithm is O(nlog 2 n). 



5.4.2 Algorithm Using Slowed Down Sorting Networks 

The above algorithm was further improved by Cole using slowed down sorting 
networks. 

A brief outline of his technique is as follows. A comparator in the sorting network 
is either active or inactive. An active comparator has its inputs determined but the 
outputs not determined. The other comparators are inactive. An active comparator 
at depth d in the network is assigned a weight l/4 rf . The intersection point of pair 
of its input, x i; , is also assigned the same weight. 

In Megiddo's algorithm we exhausted all the x,/s of each iteration before going 
to the next iteration. Here we simultaneously deal will all the x l; 's of the active 
comparators with their respective weights. Accordingly, we evaluate F at weighted 
median of x t /s in the place of median, in the search of interval that contains 
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x*. Thus, we resolve comparisons in either the left half, (-00, z m ], or the right 
half, [x m ,oo) of the weighted median point x m . These comparisons with weight 
approximately half the total weight become inactive and are replaced by active 
comparisons not more that one-fourth of the total weight as new weight. We reiterate 
these steps until we sort the lines at z*. We claim the following. 

Theorem 5.7 ([8]) At the completion of kth iteration the active weight is bounded 
by(*/ 4 ) k n/2,fork>0. 

Theorem 5.8 ([8]) For k > *>(d + l /alogn), after kth iteration there are no active 
comparators at depth d, where d > 0. 

Theorem 5.7 guaruntees that in every iteration we drop a finite fraction of total 
weight whereas Theorem 5.8 guarantees that the algorithm terminates in <9(log n) 
steps. Therefore, after O(logn) steps we know the median of / t 's at z* and hence 
x*. The running time of this algorithm is O(nlogn). 



5.4.3 Application of Prune and Search Technique 

Now we apply prune and search technique to this problem. Let us refer to Algo- 
rithm 5 presented in the previous section. The input of this algorithm is set of lines 
in the plane. We do the following: 

We prune the lines which are not the fc-th largest lines (initially k = [n/2j) 
at x* and which lie above or below the trapezium given by Theorem 5.4. 

When some of the lines are pruned we update value of k. 

The oracle H is constructed by evaluating F(p) at x = p and determining 
position of x* with respect to p with the help of sign of F(p). 

This algorithm runs in linear time which is optimal. 
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5.5 Linear Programming 

The next problem that we solve is the linear programming problem in fixed dimen- 
sions. The linear programming problem is as follows. 

Minimise z 
Subject to 

{,*) < *, l<t<n 



Each constraint defines a halfspace in a d+ 1 dimensional space. We try to compute 
the solution point x*. We again apply the Algorithm 5 for computing the fc-th 
largest hyperplane at x*. 

We construct the oracle ft for this problem as follows. We recursively compute 
the optimum value constrained on the query hyperplane H. We can compute the 
relative position of the optimal point with respect to H using the gradient of largest 
constraint at the constrained optimal point on H. The base case in the recursion is 
a point P corresponding to 0-th dimension where the solution is trivial. 

We sort the constraints and choose the largest constraint at x*. Whenever a 
comparison is resolved the smaller hyperplane is pruned. Since the highest constraint 
is to be computed, Theorem 5.5 is applicable. This algorithm runs in linear time 
though it is not as efficient as that of Megiddo because of the large constants involved 
in the AKS sorting networks. 



5.6 Ham-sandwich Cuts 

Let A,-, where 1 < i < rf, be d finite sets of points in ^. The ham-sandwich cut of 
these is a hyperplane such that it bisects each of these sets. To compute this cut 
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we dualise the problem as usual. In the dual space the sets A t transform to d sets 
of hyperplanes. We define this transformation as follows. 



Primal Space y (m,f) -f c (x, y) 

Dual Space (-m, c) Y = (x, X) + y 



It should be noted that the above transformation conserves the incidence relation- 
ships. 

A ham-sandwich cut transforms to the point, say h, such that in each set there are 
exactly half the total number of hyperplanes that are above h. Consider the problem 
of computing h. If the median hyperplanes of each set at the point h are known 
then h can be computed by computing the intersection point of these. Hence, we do 
the following in the Algorithm 5. The elements thai are to be sorted are the above 
hyperplanes and the ordering of the hyperplanes is to be done at A. A comparison 
in the sorting network would produce a "ridge". We need to localise h with respect 
to these "ridges". The comparisons corresponding to those ridges for which h is 
localised get resolved. 

In the oracle we need to know which side of a given hyperplane H does the point 
h lie. For this we use the method of Lo, Matousek and Steiger[17j. This provides 
us a method to determine if h lies inside a given infinite prism. This oracle n first 
computes the median level of one of the sets on each vertical side of prism and 
then counts the total number of intersections below this level of the median levels of 
other d - 1 sets. If this number is odd then h lies inside the prism. The complexly 
of this algorithm is bounded by the complexity of constructing a *Mevd in the 
d - 1 dimensional plane. So we get the same time bounds as in [17]. The hnearly 
separated case in E* can also be solved similarly in linear time. 
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5.7 Centrepoints in d- dimensions 

In this section we discuss the problem of computing centrepoints of sets of points 
in J-dimensional space. Suppose $ is a set of n points in 9^. The centre of S is 
union of all points p such that any closed halfspace containing p contains at least 
\nl(d+ 1)"| points. We obtain a weaker type of centre, approximate or e-centre if 
we decrease the lower bound in the above definition of centre to fn(l e)/(d + 1)] , 
where < e < 1. 

We first discuss the method of computing an approximate or e-centre of set S, where 
0<< 1. 

To compute an e-centrepoint we dualise the centrepoint problem. However, as 
opposed to ham-sandwich cut algorithm in the previous section, here we need to 
consider both the primal space as well as the dual space, in the algorithm. In 
particular, we search the centre in the dual space and prune the points in the primal 
space. As before we employ the technique by applying Algorithm 5. 

A centrepoint in the dual space will be a hyperplane (a centre-hyperplane) such 
that in every vertical line there would be at least \n/(d -f 1)] intersections of given 
hyperplanes above and below the intersection of centre-hyperplane. The dual of 
centre is union of such hyperplanes. The e-centre-hyperplane is defined similarly. 

The elements to be sorted are the dual hyperplanes of points in S. We apply sorting 
on these hyperplanes. Since We do not execute the AKS sorting network to its 
completion, we can be vague about the sorting order and so we do not specify the 
exact point at which the sorting is taking place. We apply the algorithm just to 
partition the space such that each partition contains only small number of related 
unresolved comparisons. In the absence of localisation the oracle 17 is not needed 
here. 

Whenever a hyperplane is queried in the Algorithm 5, we process the next steps 
allowing all the possible outcomes of the constrained search of the oracle ft. In this 
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way whole of the region is partitioned into several sets each containing only small 
number of intersection with "ridges" related to unresolved comparisons. We iterate 
until the maximum of these numbers drops below e* - *Cin* . Then we apply the 
prismoid method in each partition to compute 4e-centre. It can be seen that the 
region determined by \n/(d + 1)] lowest and highest levels is an 4e-centre. This 
is the consequence of Theorem 5.4. Clearly, each face of 4e-centre contains exactly 
\n/(d+l)] points. 

Since we are not searching for a solution point, we do not prune points as usual. 
We provide another method of pruning below. In the earlier chapter we saw how 
pruning can be done in the plane for collection of three and four points. These 
points were respectively deleted and substituted by the Radon's point. In fact, it 
can be shown that (d+ l)r-f 1 points in ^, for any r, can be suitably substituted by 
their Tverberg's point (a generalisation of Radon's point) when certain conditions 
are met. Unfortunately, this also can not be applied in our case. To prune the 
points, we conjecture the following[33]. 

Conjecture 5.1 If there are (d -f l)r + q, where < q < r points in E d such that 
each r-set is contained in some halfspace containing less than [n/(d + 1)J points of 
S then these (d + l)r -f q points can be substituted by q points. 

If this conjecture is proved then we will be able to discard at least a fraction of 
points in each iteration and a centrepoint in d dimensions can be found in linear 
time. 



Chapter 6 



Conclusions 



An optimal algorithm has been described in this thesis for computing the intersec- 
tion radius of a given set of planar objects containing points, straight lines, rays, 
line segments, wedges, half planes and planes. We have used a new approach of 
allowing addition of objects in the pruning process while ensuring a simultaneous 
net reduction of the input set. This provides one more example of the power and 
versatility of the prune and search paradigm. This algorithm has applications in 
many practical stabbing radii problems. 

We have also presented an optimal algorithm for computing a centrepoint for a 
planar finite set of points. This settles a long standing problem in computational 
geometry. It may be worthwhile to point out that an approximate centrepoint 
suffices in many applications but this exact algorithm, because of its efficiency, can 
very well replace the existing approximate centrepoint algorithms in two dimensions. 

A new technique has been presented to solve geometric optimisation problems using 
prune and search paradigm. We have been able to do this by the synthesis of para- 
metric searching and prune and search techniques. We have applied this technique to 
compute ham-sandwich cuts and centrepoints in higher dimensions. This technique 
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has also been applied to the well known problem of linear programming. The results 
obtained match the best results extant for these problems. 



6.1 Further Research Problems 

It would be worth investigating if there exist similar algorithms as presented in this 
thesis to compute intersection radius of objects in higher dimensions. The weighted 
version of this problem also remains unsolved. The intersection radius for a set 
of simple polygons cannot be computed using the same algorithm because of the 
inherent non-convexity of the distance function. 

Just as there are centrepoints, we can have other centre objects like centreline, 
centrecircle etc. The computation of these remain an open problem. The weighted 
version of the centrepoint problem also remains open. The other problems for a 
finite set of simple polygons in the plane to look at are area-centrepoints, perimeter- 
centrepoints and the like. The weighted centrepoint can be addressed as a special 
case of these problems. 

In this thesis, a centrepoint of an explicit set of points is computed in linear time. 
However, the centrepoint can not be computed within the same time bounds if the 
point set is implicit, such as the set of intersection points of straight line defined by 
any two points of a set of points in ft 3 with a fixed plane. 

In the technique of parametric prune and search the open area of research is to 
provide efficient algorithms for pruning. The pruning technique described in this 
thesis has a doubly exponent factor of d in its time complexity which renders 
the algorithm impractical for higher dimensions. AKS sorting networks have also 
large constants of proportionality in their complexities. In practice, to solve the 
problem discussed in this thesis faster, randomised algorithms can be used. The 
other important question is how this technique can be applied to other parametric 
searching algorithms as well. It seems that the method presented in the thesis is a 
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particular case of Mataousek's result [19]. It remains to be seen if this technique can 
be applied to non-linear programming problems as well. 

Furthermore, AKS sorting networks have built-in comparators that take two inputs 
each. The generalisation of comparison for points in higher dimensional space is their 
rotational order. So, if a network is designed which has such rotational ordering 
elements in place of comparators then it can be applied to many more problems 
using a similar technique. 
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